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ABSTRACT 


The flowfield about a Vertically-Launched Surface-to-Air Missile model 
at an angle of attack of 50° and a Reynolds number of 1.1x105 was investi- 
gated in a low speed wind tunnel at the Naval Postgraduate School. The goal 
of this thesis is to determine the location and intensity of the asymmetric 
vortices in the wake of the VLSAM model and to display these vortices by 
velocity mapping and pressure contours. The two model configurations tested 
were for a cruciform missile with wings and tails; one at 0” roll angle ("plus" 
aspect) and the other at a 45? roll angle ('cross" aspect). Two flowfield 
conditions were treated: the nominal ambient wind tunnel condition and a 
condition with a grid-generated turbulence of length scale 1.08 inches and 
1.88% turbulence intensity. The turbulence length scale is 61.7% of the model 
diameter and 4.7% of the model length. The following conclusions were 
reached: 1) An increase in turbulence intensity tended to reduce the strength 
of the asymmetric nose-generated vortices; 2) the two asymmetric vortices 
remained in approximately the same position for an increase in turbulence; 3) 
"cross. aspect vortices were more diffused, slightly larger and centered 
further away from the model surface than those of the "plus" aspect body 
configuration, which correlates with the differences in induced side forces for 
these configurations observed by Rabang; 4) the top vortex of the two 
asymmetric vortices was closer to the model surface and appeared to be 
stronger for both configurations; and 5) the addition of wings and tails did not 


greatly alter the vortex pattern around the nose of the missile model. 
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I. INTRODUCTION 


A. BACKGROUND 

The development of the Vertical Launch Surface-to-Air Missile (VLSAM) 
has provided greater weapon system reliability, availability and flexibility 
over its predecessors. Earlier systems, such as trainable launchers and box 
launchers, were rather cumbersome and required excessive deck space, thus 
limiting the number of missiles per ship. Additionally, these systems were 
slow to reload since trainable launchers required cycling of rounds to get the 
desired one in position for launch, while box launchers were reloaded by 
hand. The vertical launcher provides for quick access to any round without 
cycling, rapid reloading and, due to its design and the fact that the missile 
blast is kept in a concentrated area near the launcher, it requires less deck 
space. [Refs. 1, 2 and 3] 

The VLSAM 's trajectory allows it to point to its target after launch and 
subsequently guide itself to the correct heading. Its aerodynamic characteris- 
tics may significantly change as it operates from subsonic to supersonic 
speeds during the launch, midcourse and terminal phases of its flight. [Ref. 2] 
The varying flight control requirements and the airflow about the missile 
during these phases provide for different flight regimes. 

When launched, the VLSAM enters an open ocean environment and is 
subject to potentially significant cross-winds, the result of which 1s a missile 
flying at relatively low velocities at a high angle of attack. [Refs. 1 and 2] In 
particular, an example of this low velocity/high angle of attack condition is 


the Standard Missile 2-Block 4, or AEGIS Extended Range version, which 


has an added sustainer section and leaves its launcher at a much lower speed 
than the unboosted SM-2 version. À missile flying at relatively low velocities 
at a high angle of attack may bring with it the formation of asymmetric 
vortices about the missile nose and afterbody which can generate unpre- 
dictable and undesired side forces, thereby affecting overall missile flight sta- 
bility. The airflow around a ship's superstructure, the ocean surface condi- 
tions and the atmospheric boundary layer conditions during launch may also 
provide turbulent flow over the missile. [Ref. 2] The advent of the VLSAM 
and the desire to have highly maneuverable supersonic missiles have 
increased the need for further studies in high angle of attack research under 
various flowfield conditions. 

This thesis continues experiments and research conducted to date at the 
Naval Postgraduate School (NPS) to understand the effects of turbulence on 
the VLSAM aerodynamic characteristics. Previously, Roane [Ref. 1] developed 
a system model to measure flowfield turbulence through the use of four dif- 
ferent grids which generate varying turbulence levels in the NPS low speed 
wind tunnel. Rabang [Ref. 2] studied the effects of this turbulence on the 
asymmetric vortex forces on the missile model. Lung [Ref. 3] determined the 
location and intensity of the asymmetric vortices in the wake of the model by 
experimental flowfield measurements about a body-only missile configuration 
with and without freestream turbulence. Similarly, the goal of this study is to 
determine vortex locations and intensities, through a series of wind tunnel 
experiments, for two missile configurations, both with and without turbu- 


lence. The configurations considered are for a cruciform missile with wings 


and tails; one at 0° roll angle ("plus" aspect) and the other at a 45? role angle 


("cross" aspect). 


B. HIGH ANGLE OF ATTACK AERODYNAMICS 

In high angle of attack aerodynamics, flow separation from the body, 
wing and tail surfaces is important due to the strong normal and side forces 
which may be generated. Major factors which influence flow separation are 
nose shape, angle of attack, crossflow Reynolds number and nose fineness 
ratio. Other factors include roll angle and roll rate, free stream turbulence, 
surface roughness, acoustic environment and VLSAM model vibrations 
(Ref. 4). 

1. Aerodynamic Regimes 

As the angle of attack, =, of a slender body of revolution ranges 

from 0° to 90°, there are at least four distinct airflow regimes through which 


the missile body transitions. [Ref. 4] 


(1) Regime I (0°< œ < 5°): At very low angles of attack, the axial flow domi- 


nates and there is no discernable boundary layer separation (flow is 
attached). 


(2) Regime II (5?« œ < 20°): At intermediate angles of attack, boundary 
layer separation occurs on the lee side of the body. This becomes a free 
shear layer which rolls up into a symmetric vortex pair that is steady 
with time. No side forces or yawing moments are induced. 


(3) Regime III (20?« « « 60?): At high angles of attack, crossflow effects 
start to dominate and vortices are now shed asymmetrically. These 
vortices induce side forces (out-of-plane forces) and yawing moments. 
The more asymmetric the vortex, the greater the side force magnitude. 
There are some flow instabilities toward the higher end of this angle of 
attack range. 


(4) Regime IV (60°< = < 90°): At very high angles of attack, the crossflow 
completely dominates and flow separation is unsteady. The Reynolds 
number, Mach number and geometry determine whether the boundary 
layer is shed as a von Karman vortex street or a random wake-like 


flow. [Refs. 1, 2 and 3] 


Figure 1 shows sketches of all four vortex-generation regimes. 
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Figure 1. Airflow Regimes [Ref. 2] 


2. Asymmetric Vortex Theory 

The principal cause of the formation of asymmetric vortices is still 
not completely understood and may be attributed to many factors. One idea is 
that boundary-layer-induced asymmetry in the location of flow separation 
causes the vortex flowfield to become asymmetric. These boundary layer 
asymmetries may result from transition and separation differences on each 
side of the missile body. Another proposition is that a hydrodynamic (inviscid) 
instability in the pair of initially symmetric vortices causes the asymmetry. 
[Refs. 5, 6 and 7] These vortices, which increase in strength with angle of 
attack, interact with the surrounding potential flowfield to provide the 
asymmetric configuration. A vortex-switching phenomenon has also been 
observed in which the vortex pattern rapidly switches from an almost sym- 
metric to a highly asymmetric configuration, which may possibly relate to a 
second inviscid solution in the leeward flowfield. [Refs. 8 and 9] At any rate, 
even though their major cause has not been determined, the behavior of 
asymmetric vortices has been well documented for a large number of models 

and shapes. 
Nose-generated asymmetric vortices appear in the flowfield around an 
ogive-nosed, slender, cylindrical body. The vortex formation occurs along the 
entire body length and induces a significant side force on the body. Figure 2 


shows this vortex flow along the length of an unyawed, slender nose cylinder. 
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Figure 2. Vortex Flow About a Slender Nose Cylinder [Ref. 1] 


3. Two-Dimensional Crossflow 

Airflow over a slender body can be divided into normal and axial 
components. Axial flow is along the missile body length while crossflow is 
essentially a two-dimensional flow normal to a cylinder. The crossflow 
analogy provides information for cylinder lift and drag which act in the 
crossflow direction. Depending on the type of flow separation on either side of 
the cylinder, side forces (at right angles to the crossflow) may exist. [Ref. 2] 

Boundary layer transition and separation mechanisms may provide 
an explanation for flow separation and subsequent asymmetric vortex genera- 


tion. The primary factor which influences the separation location of the 


boundary layer is the crossflow Reynolds number. Other factors include sur- 
face roughness and turbulence. 

Flow around a cylinder in incompressible flow can be classified into 
four regions, represented by differing flow separation and drag behavior, as 
depicted in Figure 3. [Refs. 10, 11 and 2] In the subcritical range, the bound- 


ary layer is laminar and flow separation occurs close to the lateral meridian 
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Figure 3. Two-Dimensional Crossflow About a Cylinder [Ref. 2] 


where the angle ($6) from the crossflow direction varies from 80° to 90°. 
[Ref. 12] In the critical Reynolds number region, a drag bucket is produced. 
The laminar boundary layer separates from the body at 6 = 90° and reat- 
taches as a turbulent boundary layer which is more energetic. Separation is 
delayed to o = 140°, resulting in a reduction of the drag. [Ref. 3] A laminar 
bubble is formed between the laminar separation and the turbulent reat- 
tachment. At this point, the flow separation may easily fluctuate from critical 
to subcritical for small changes in Reynolds number. From Figure 4 [Ref. 10], 
when critical separation exists on one side of the body while subcritical sepa- 
ration is on the other, a large difference in 6 is possible. Therefore, vortices 
will be at maximum asymmetry and the side force at the highest magnitude. 
Since maximum side force occurs at the critical Reynolds number, it 1s a 
noteworthy parameter. [Ref. 13] 

As the Reynolds number is further increased, turbulent separation 
moves forward to ὁ « 140° and the laminar bubble no longer exists. At 
$ > 140°, asymmetric vortices are ineffective at producing a significant side 
force, thus the sudden decrease in magnitude (Cy/Cy) as shown in Figure 4. 
For supercritical and transcritical Reynolds numbers, the laminar transition 
point moves towards 6 = 0° and turbulent separation occurs at 6 = 100°. The 
asymmetric transcritical separation point moves towards the lateral merid- 
ian, where the vortices once again produce a significant side force. [Refs. 5 
and 2] The Reynolds number provides the greatest influence on the normal 
force and drag characteristics, particularly within the critical range where 


the maximum normalized side force and maximum vortex asymmetry occur. 


[Ref. 12] 
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Figure 4. Side Force to Normal Force Ratio [Ref. 5] 


Another study by Lamont [Ref. 14] describes a different effect of 
Reynolds number on the maximum side force as illustrated by Figure 5, 
where the side force at an angle of attack of 55° is plotted. The maximum side 
force falls from a high value at laminar separation to almost zero in the 


middle of the transition region before climbing again to higher levels at fully 
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Figure 5. Effect of Reynold's Number of Maximum 
Side Force at « = 55° [Ref. 14] 


turbulent separation. Thus, there appears to be two different mechanisms for 
producing asymmetric flow and, hence, a side force on an ogive-cylinder. One 
mechanism operates in both the laminar and the fully turbulent separation 
regimes, in which the side force results from asymmetric vortex patterns in 
the wake of a body. The other mechanism occurs only in the transitional 
separation regime. Here, the Reynolds number at which the near-zero side 
forces were recorded, is the same range of Reynolds number in which the 
minimum drag coefficient on a 2-D cylinder occurs and in which no coherent 


vortex shedding can be detected. 
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4. Three-Dimensional Vortices 

The missile nose geometry is an important factor in vortex genera- 
tion and disposition since vortices shed at the nose tend to dominate other 
vortices along the body length. [Ref. 15, 8 and 16] 

Nose-generated vortices are sensitive to the nose roll angle due to 
surface imperfections and nose geometric deviations [Ref. 17]. Rabang varied 
the roll angle in 45° increments and investigated the resulting side force coef- 
ficients, shown in Figure 6. The vortex system generated by the nose domi- 
nates afterbody vortices for body configurations with and without wings 


regardless of the turbulence conditions [Refs. 2, 18, and 19]. 
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Figure 6. Side Force Variations With Nose Roll Angle [Ref. 2] 
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Nose geometry may be pointed or blunt for cones and ogives. For 
pointed noses, angle of attack for the onset of asymmetric vortices (=y) is a 
function of the semi-vertex angle (64). Asymmetric vortices start at the nose 
and are rapidly shed, yielding unsteady side forces. [Ref. 3] At all Mach 


numbers, asymmetric vortex shedding starts when «ay is greater than 04. 


For a conical nose: 


04 7 e Av/ 2 (approximation) (1) 


For a tangent ogive nose: 


d EE 


- tan 
A 2 
1 4-0.25 
Ν (2) 
where lw is the nose length and d is the base diameter, or 
6,4 = ln/d (approximation) (3) 


Nose fineness ratio also affects the asymmetric vortex induced side 
forces in that as this ratio increases, the side force also increases. With an 
increasing ratio, both the nose apex angle (04) and the angle of attack for the 
onset of asymmetric vortices (e«Avy) will decrease, making the missile more 
susceptible to induced side forces at lower angles of attack. [Refs. 15 and 20] 
Decreasing nose fineness ratio has been found to be more beneficial in reduc- 
ing side forces than blunting the nose [Ref. 15]. Side force decreases with an 
increasing Mach number to nearly zero at supersonic speeds [Refs. 18 and 


21). 
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C. TURBULENCE 

Turbulence denotes the presence of random, short duration variations in 
a flowfield with a given mean velocity. When calculating turbulence effects on 
a body in the flowfield, a comparison between the scale of the body and that of 
the turbulence must be made. The energy in the turbulence flowfield should 
also be considered. [Ref. 1] 

Turbulence intensity, T,, is the measure of the relative magnitude of 
velocity fluctuations in the flowfield. For a horizontal flowfield or crosswind, 
it is the ratio of the root-mean-square (rms) velocity fluctuation, u', to the 


mean velocity component in the flowfield, Uee. 
Tu 2 u/Uee (4) 


Turbulence length scales describe the time-averaged measure of the size 
of the constantly changing fluid disturbance eddies. Àn increase in the spatial 
length of the turbulence corresponds to an increase in the time the body is 
exposed to the fluctuation. Large and small scale turbulence length scales are 
both found in a flowfield. 

From a single source, the cascade effect produces turbulence eddies of 
different length scales. This "cascade" effect is caused by a strain in one direc- 
tion (x, y, z plane) which affects the orthogonal components due to the conser- 
vation of angular momentum. For example, an increase in the x and y velocity 
components of a vortex rotating in the x-y plane will have an effect on the ve- 
locity and length scales of the y-z components. [Ref. 1] Cascading continues 
until the smaller eddies disappear due to the viscosity. As turbulence 


decreases, the energy transfer decreases and the individual intensities of 
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each eddy will decrease at a faster rate. [Ref. 22] Thus, the larger scale turbu- 
lence predominates. 

The length scale to body size ratio may determine the manner in which 
turbulence affects the VLSAM flowfield. It may be compared to missile 
length, L,:L), or missile diameter L,:Lg. [Refs. 1 and 2] 

For length scales much greater than the body, L,>>L}, the effect is like a 
steady-state flowfield, where deviations in speed and direction would be of 
long duration. The flowfield effects on vortex development are mainly domi- 
nated by the same factors and conditions as for a two-dimensional cylinder. 

In contrast, unwanted rolling, pitching and yawing motion of the body is 
primarily caused when the turbulence length scale is comparable to the body 
length, Ly=L). [Ref. 23] The flowfield is distinctly non-steady for this case. 

When the length scale is of a dimension much smaller than the body, 
most significantly, when it is smaller than the missile diameter, Ly<<lLg, it 
has a magnitude comparable to the boundary layer thickness on the missile 
surface. Thus, boundary layer development and flow separation over the body 
may be affected by the presence of smal] scale turbulence. An increase in tur- 
bulence intensity with a length scale on the order of the boundary-layer scale 
tends to reduce the magnitude of induced side forces. [Refs. 2, 24, and 25] A 
goal of the current investigation is to determine the effect of vortex-scale 


turbulence on the asymmetric vortices and resulting induced side forces. 


D. LIFTING SURFACE EFFECTS 
The complete vortex structure of the missile is a net result of the indi- 
vidual contributions from the body, wings, strakes and tails. In general, mis- 


siles use low aspect ratio wings (when compared with aircraft). Since some 
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missiles have wing spans that approach body diameter, it is important to 
consider the joint effects from a wing-body combination. Nose vortices dictate 
flow behavior over a missile body at high angles of attack and, consequently, 
these vortices may also be felt by the wings. Nose and body vortices move 
away from a missile body without wings but, when wings are added, they 
move closer to the body. This result is comparable to increasing the effective 
angle of attack causing unsteady asymmetric vortices. For wings with low 
aspect ratio, a major portion of the lift produced by the wing will be a result of 
vortex lift. The net effect of the wing-body combination appears to be a reduc- 
tion in the effective angle of attack for the onset of asymmetric vortices and 
side forces. [Refs. 2, 3, and 26] 

Vortex lift effects are improved by incorporating strakes with low aspect 
ratio wings. The strakes produce additional strong vortices. Some researchers 
have found that placing long strakes on a missile would induce interference 
with the crossflow component around the body, thus decreasing the effect of 
the forces and moments generated by asymmetric body vortices. [Ref. 6] 
Rabang has shown that the addition of typical VLSAM wings and strakes 
tend to preserve the induced side force for all levels of turbulence intensities 
and length scales. [Ref. 2] 

The addition of tails has a minor influence on forebody flowfields and 
maximum side forces, in particular, at low angles of attack. At higher angles 
of attack, nose and wing vortices may have a slight effect on the tailflow, 
depending on the wing placement and afterbody length between the wings 


and tails. 
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E. VLSAM LAUNCH ENVIRONMENT 
1. Marine Environment 

Turbulence conditions which exist within the atmospheric boundary 
layer (ABL) may significantly impact the VLSAM. The ABL is the lowest por- 
tion of the atmosphere and is formed by its interaction with the surface over 
which it flows. Turbulence in this layer is the result of the transfer of heat, 
momentum and mass. 

The surface layer, the lowest segment of the ABL, can vary in 
height from 5 to 200 meters but is typically on the order of 50 meters. It is 
also characterized by mechanically produced, small-scale turbulence resulting 
from surface roughness or friction from waves on the ocean surface. This 
small-scale turbulence is larger than the missile length. This region is 
described by variations in wind speed, nearly vertical heat and mass fluxes, 
and other meteorological fluctuations with height. [Refs. 27 and 23] Further- 
more, the majority of the flow in the surface layer itself can be considered 
horizontally homogeneous and two-dimensional (Ref. 23). 

A measure of the roughness of the surface is called the roughness 
length, Zo, which is a function of the mean wind velocity at various heights 
above the surface. By combining the roughness length with the elevation and 
wind speed, both the turbulence intensity and length scale can be empirically 
determined. [Refs. 2, 3 and 28] Typical open ocean surface roughness lengths 
are in the range of 0.001<Z,<0.01, with Z, in meters. For a 10 meter elevation 
and a mean wind speed of 25 m/sec in a neutral atmosphere, turbulence 
intensities may range from 13 to 17 percent. [Ref. 28] The longitudinal 
turbulence length scale would then range from 262<L,<295 feet (80<Ly< 90 
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meters). Therefore, for a typical missile with a 1.1 foot diameter, the 
turbulence length scale to missile diameter length scale ratio is Ly:Lg=280:1. 
[Ref. 28] This represents a length scale very much larger than a conventional 
missile length and, therefore, would have little effect on its boundary layer 
development. However, the cascade effect from large scale turbulence and 
from crosswind interaction with a ship's superstructure allows length scales, 
initially much larger (85 meters) than the dimension of a missile, to decrease 
(cascade) to scales where they could affect the missile boundary layer 
development and the development of vortices from the missile nose. The 
actual amount of such small-scale turbulence present in the marine 
atmosphere is largely unknown, however. [Ref. 23] 
2. Launch and Crosswind Velocities 

A typical VLSAM at launch (vertical velocity = 164 ft/sec) is still 
well within the surface layer environment and is subject to both crosswind 
and turbulence effects. [Ref. 29] Crosswind velocity depends on both the 
ambient wind speed and the speed of the launch platform. A ship speed of 20 
knots with a head wind of 20m/sec, combined with the VLSAM launch 
velocity, results in an effective angle of attack of 31° at 191 ft/sec, which 
places the missile in the asymmetric vortex region (Regime III) almost imme- 
diately after clearing the exit plane of the launcher, 0.2 seconds after launch. 
[Refs. 1 and 2] 

When the missile first leaves its launcher it will experience an even 
greater effective angle of attack due to its slower velocity. A ship's hull and 
superstructure can dictate changing flow fields and turbulence at this initial 


launch altitude. The ship airwake may increase crosswind velocities and 
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cause significant crosswind gradients in the flowfield, thus increasing turbu- 
lence intensities while decreasing turbulence length scales. [Ref. 2] 

Later in its launch profile, when the VLSAM pitches over towards a 
target, it may reach effective angles of attack of up to 50?. [Refs. 2 and 30] 
Thus, during the launch phase there exists a definite possibility of asymmet- 
ric vortex induced side forces on a VLSAM. 

3. Additional Launch Considerations 

During launch, a missile is influenced by many factors which are 
dependent on missile design, ship's orientation, launcher mechanics and 
missile flight control. Shipboard roll, pitch and yaw are directly transmitted 
to the launch platform and must be taken into account. Inherent factors, such 
as the plume (jet) effect of the missile's engine and blast effects of the vented 
exhaust gases, can also affect VLSAM aerodynamics. Exhaust gases can 
impinge directly on the missile surfaces or they can impact the flowfield in 
which the missile 1s launched, especially if the gases are vented upward into 
the vicinity of the accelerated missile. [Ref. 1] The manner in which control 
systems respond to missile orientation changes is another factor. Should the 
missile change its flight attitude, the flowfield around it will also be altered. 
Obviously, the considerations discussed in this section are not all-encompass- 
ing. There are many other factors that affect missile flight behavior during 


launch, but they are beyond the scope of this thesis and will not be included. 


18 


II. EXPERIMENT AND PROCEDURES 


A. PURPOSE 

In this study, the location and intensity of asymmetric vortices in the 
wake of the VLSAM model were determined for varying levels of turbulence. 
This turbulence was generated by the placement of a series of grids in the 
wind tunnel. The vortices were displayed by velocity mapping and pressure 
contours. To accomplish this, wind tunnel flowfield pressure measurements 
were taken for a specified survey grid using a scanivalve/probe data acquisi- 
tion system. 

Figure 7 [Ref. 3] shows the planar survey grid, the x-y plane, which was 
perpendicular to the freestream velocity and located 10.5 inches downstream 
from the missile model's nose. The model body was centered on the y dimen- 
sion. For the actual data acquisition runs, 23 points were measured in the y 
direction. There were 11 points above and below the model centerline, with 
point 12 directly at the centerline. Along the x axis, 18 points were measured. 
the x-y dimension for the experiment was 3x5.5 inches, with a step distance 
(increment) of 0.25 inch (22x0.25=5.5 and 12x0.25=3 inches). This dimension 
covered the main portion of asymmetric vortices. 

Pressure measurements were obtained by the 5-hole probe throughout 
the survey plane. The data from the pressure probe was reduced through the 
use of computer programs to obtain isobars of total pressure coefficient and 


static pressure coefficient, and to map the crossflow velocity vectors. These 
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Figure 7. The Planar Survey Grid [Ref. 3] 
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results were correlated with the force measurements of Rabang and with the 
previous experiments by Lung to provide a greater understanding of the 
vortex flowfield. The following sections further discuss the equipment and 


software used, and the experimental procedures followed. 


B. APPARATUS 
Information about the construction, specifications and configurations of 
the major pieces of equipment used in this study is described in this section. 
1. Wind Tunnel 

The low-speed, single return, horizontal-flow tunnel located in 
Halligan Hall at NPS was utilized. (Figure 8, [Ref. 31]) . It is powered by a 
100 horsepower electric motor coupled to a three-blade variable-pitch fan via 
a 4-speed Dodge truck transmission. Aft of the fan blades are a set of stator 
blades which help straighten flow. Two fine wire mesh screens located up- 
stream of the settling chamber plus turning vanes at all four corners reduce 
turbulence. A heavy wire screen behind the test section prevents foreign 
object damage to the fan blades. [Refs. 1, 2 and 31] The tunnel is 64 feet long 
and ranges from 21.5 to 25.5 feet wide. 

The wind tunnel test section measures 45 inches by 32 inches. The 
walls diverge slightly to prevent reduction in freestream pressure due to 
boundary layer growth. The settling chamber area to test section area 
contraction ratio is approximately 10:1. Corner fillets, which house the light- 
ing, reduce the section area from 10 ft? to 9.88 ft?. Fillets are found at wall 
intersections throughout the tunnel to help reduce boundary layer effects. 


[Ref. 31] 
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Figure 8. Naval Postgraduate School Wind Tunnel [Ref. 31] 


A reflection plane installed in the test section reduces the available 
height to 28 inches. A flush-mounted turntable allows for changes in model 
pitch angle or angle of attack via a remotely controlled electric motor beneath 
the tunnel. Since the test section operates at atmospheric pressure, breather 
slots are installed around the tunnel perimeter to replenish air lost through 
leaks and to ensure a uniform test section pressure. The tunnel was designed 


to provide test section velocities of up to 290 ft/sec. [Ref. 31] 
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Wind tunnel temperature is measured by a dial thermometer 
extending into the settling chamber. Dynamic pressure is measured by the 
static pressure difference between the test section and the settling chamber 
using a water filled manometer. The static pressure is measured by four pres- 
sure taps located upstream from the test section to avoid interference from 
the model. These taps are connected via a common manifold prior to feeding 
into the manometer. Pressure differences are measured in centimeters of 
water. Equation (5) is used to convert to the actual wind tunnel velocity. 


[Ref. 3] 


(22.0475) (P_H,0) ] 
υ -ἰ---------- (5) 


(K) (p) 


where: 


Um = measured velocity (ft/sec) 
2.0475 = conversion factor 
Pem H20 = manometer reading 


K = calibration factor (for specific grid) 
p = air density (lb/ft?) 


2. Turbulence-Generating Grids 
Four grids are used to create turbulence of varying intensities and 
length scales. Each is mounted in a wooden frame and placed 73 inches for- 
ward of the pivot axis of the model support system (see Figure 9). Three of the 
grids are constructed from wood and the fourth is made of wire. Grid specifi- 


cations are listed in Table 1 and are shown in Figure 10. They are 
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Figure 9. Planview of VLSAM Model With Pressure Probe and Grid In the 
Test Section of the Wind Tunnel (not drawn to scale) [Ref. 2] 
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TABLE 1. GRID SPECIFICATIONS [REF. 2] 


Grid Mesh Width (in.) | Bar Diameter (in.) | Nfesh/Diameter Material 


Four 1.00 0.0625 16 Wire 





Figure 10. Square-Mesh Turbulence-Generating Grid [Ref. 3] 
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square-mesh square-bar biplanar grids which generate nearly isotropic homo- 
geneous turbulence. (Ref. 29] Roane measured turbulence intensities and 
estimated length scales, shown in Figures 11 and 12. (Ref. 1] The grid turbu- 
lence parameters taken by Roane are summarized in Table 2. Grid turbulence 
effects, with respect to changing length scales at constant intensity or con- 
stant length scales with changing intensities, can not be investigated with the 


present grid geometries. [Ref. 3] Figure 13 shows photographs of the grids. 
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Figure 11. Grid Turbulence Intensities (dashed line indicates 
model pivot axis) [Ref. 2] 
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Figure 12. Grid Turbulence Length Scales [Ref. 2] 


TABLE 2. GRID TURBULENCE PARAMETERS 
(AT MODEL PIVOT AXIS) [REF. 1] 









Intensity Length Scale Turbulence/ Dynamic Pressure 


Grid (percent) (in.) Model Dia. (Ib/ft?) 


None 0.23 15.85 
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Figure 13. Turbulence-Generating Grids 
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3. VLSAM Model and Support Equipment 

The model vvas designed to represent a current cruciform tail-con- 
trol missile with very low aspect ratio wings (long dorsal fins). It was con- 
structed from 6061 and 2024 aluminum alloy by NPS personnel. [Ref. 1] 

The hollow cylinder body section contains locating pin attachment 
points for the balance, sleeve, wings and tails. The machined sleeve provides 
a close tolerance fit between the balance gage and the interior of the model. 
[Ref. 3] Both body roll angle and nose roll angle may be varied in 45° incre- 
ments. The wings with strakes and the tail control fins are rigidly connected 
to the model body by countersunk screws. Figure 14 depicts the dimensions 
and specifications of the VLSAM model. [Ref. 2] The model's surface is 
polished and free of protruberances. 

The model support is ngidly fixed in the test section by the reflec- 
tion plane turntable at the base and an aluminum reinforced clear plexiglass 
section at the top. The pivot point of this rotating support coincides with the 
approximate center of the VLSAM model. The plexiglass has three slots (7-, 8- 
and 10-inches long) cut in it, each 5/4 inches wide. These slots correspond to 
the positions of model length to diameter ratios of 3, 6 and 9; i.e., 5.25, 10.5 
and 15.75 inches from the nose. [Refs. 2 and 3] 

4.  Velmex 8300 3-D Traverser 

The Velmex 8300 is composed of a motor controller assembly and a 
traversing assembly, and uses three microcomputer-controlled stepping 
motors (one for each axis of movement). The motor controller assembly is 


capable of interpreting motor movement commands from a host computer, 
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Total length = 22.85 in. 
Base diameter = 1.75 in. 
Length/diameter ratio = 13.06 
Ogive nose length = 4.0 in. 


Ogive/diameter ratio = 2.29 


Wing span/root chord = 3.13 in./13.55 in. 
Tail span/root chord = 5.50 in./1.70 in. 
Center of pressure = 13.5 inches aft of nose tip (approx.) 
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Figure 14. Specifications of VLSAM Model [Ref. 2] 


programmable control or terminal. Software commands allow the operator to 
select motor variables such as velocity, acceleration, increment distance and 


units (motor steps or inches). [Ref. 32] 
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The stainless steel and aluminum traverser assembly (Figure 15) 


consists of three separate motor/jackscrew assemblies. 





Figure 15. Velmex 8300 Traversing Assembly 


The traverser was mounted to existing hardware on top of the tun- 
nel so as to minimize tunnel-induced vibrations. A 5-hole pressure probe, 
attached to the 8300 control drive, can be accurately and effectively moved 


through the test section. 
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5. 5-Hole Pressure Probe 

The three-dimensional 5-hole probe, Figure 16 [Ref. 33], is made of 
corrosion resistant non-magnetic stainless steel. It is 0.125 inch in diameter 
and 24 inches in total length with 22 inches of reinforcement tubing. The 
probe has five measuring holes located on its prism-shaped tip. A centrally 
located hole (P1) measures total pressure, while two lateral pressure holes 
(Po, P3) are used to determine yaw angle of flow. Pitch angle is determined by 
pressure holes (P4, P5) located above and below the total pressure hole. The 
probe is usable for speeds up to Mach 0.7. 

The speed of reading depends on the length and diameter of the 
pressure passage inside the probe, the size of the pressure tubes to the 
manometer, and the displacement volume of the manometer. [Ref. 33] For 
smaller diameter tubes, the time constant increases rapidly. For this experi- 
ment, the tube diameter was 1/4-inch O.D. and the tube lengths were three 
feet, so the time delay was about 0.15-0.26 second. [Ref. 3] 

6. Scanivalve and HP Data Acquisition System 

One 48-port scanivalve was used to measure each of the 5-hole 
probe pressures. The Hewlett-Packard (HP) data acquisition system consists 
of a combination of hardware and software that enables the IBM PC-AT com- 
puter to act as a fully automated instrumentation system. [Ref. 34] Individual 
instruments include the Relay Multiplexer, Digital Multimeter and Relay 


Actuator. 
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Figure 16. The 5-Hole Pressure Probe and Measuring Tip [Ref. 33] 


The scanivalve mechanism puts out a 7-bit binary coded decimal 
(BCD) signal that corresponds to the port (1-48) currently connected to the 
scanivalve transducer. This allows remote electronic monitoring of the port 
assembly configuration. [Ref. 3] The scanivalve consists of a transducer, 
motor drive, port assembly and solenoid controller, which actually regulates 


the scanivalve. Two commands allowed by this solenoid are STEP, which 
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moves the scanivalve to one port location, and HOME, which sends the scani- 
valve to port number 48. The Relay Actuator is used solely in controlling the 
scanivalve to STEP or HOME. [Ref. 3] 

The HP Data Acquisition System is shown in Figure 17. The scanivalve 
signal, containing probe port pressure information, passes to the Relay Multi- 
plexer which provides one common output channel for the Digital Multimeter 
(DMM). The signal is conditioned by a low pass filter prior to being measured 
by the multimeter. The DMM automatically converts input analog voltage 


signals into a digital (binary) form which can be read by the computer. 


C. EXPERIMENTAL CONDITIONS 
To facilitate data correlation, the conditions for this experiment were 
similar to the conditions of the previous studies by Lung and Rabang. 


(1) Test section reference dynamic pressures were set at 7.2 cm H20 for the 
no grid run and at 10.0 cm H90 for grid #3, which yielded a subcritical 
Reynolds number of Re=1.1x105. These reference pressures are the 
same as those used in the turbulence mapping by Roane [Ref. 1] and 
were duplicated in this study in order to ensure comparable test 
section velocities and turbulence grid length scales and intensities. 


(2) The VLSAM model nose geometry was held fixed at nose position eight, 
which Rabang showed gave the maximum side force magnitude. 
[Ref. 2] 


(3) Afterbody roll angles were as follows: 
° Body A: wings and tails at roll angle dR=0° in a '+' configuration 
* Body B: no wings or tails at roll angle ¢R=45° 


e Body C: wings and tails at roll angle óg245? in a "x" configuration. 


Figure 18 shows these three configurations. Only Bodies A and C were 
tested in this study. 
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(4) 


(5) 


BODY B 





Figure 18. VLSAM Model Body Configurations [Ref. 2] 


Model blockage factor corrections were calculated by Rabang for each 
body configuration. These factors are a function of the model angle of 


attack. For this experiment, the angle of attack was fixed at 50°, where 
the total blockage correction € equalled 0.0123. This factor was imple- 
mented in data conversion programs. 


The longitudinal position for data acquisition was at a length/body 
diameter ratio of 6, which was 10.5 inches from the nose of the missile 
model. 


(6) Wind tunnel temperatures were not allowed to vary by more than 20°F 


from the beginning to the end of a run. Wind tunnel settling chamber 
temperatures tended to rise quickly due to air friction, particularly 
when the grid was added. When temperatures were excessive, the tests 
were stopped and the air in the tunnel was circulated until it cooled 
down sufficiently before tests were continued. 
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D. SOFTWARE AND PROCEDURES 

In order to correlate results with previous data, the computer programs 
used (or developed) by Lung to acquire and reduce data were also used for 
this study. Figure 19 [Ref. 3] is a schematic flowchart of the various programs 
and their resulting data files. The following sections provide further elabora- 


tion on these programs. 
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Figure 19. Program/Data File Flowchart [Ref. 3] 
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1. PPROBE Program (Data Acquisition) 

The BASICA application program which runs the VLSAM experi- 
ment is comprised of STATEFILE, PGMSHEL and traverser programs. 
STATEFILE gives the computer configurations of the data acquisition 
instruments, while PGMSHEL informs the computer of all the functions 
available at each of these instruments. STATEFILE and PGMSHEL perform 
initialization chores and allow communication between the HP instruments 
and the IBM computer. [Ref. 34] 

The traverser program is the actual application code which allows 
the operator to precisely control pressure probe movement by either manual 
or computer-controlled input. The traverser program was written by 
Kindelspire [Ref. 35] in the Advanced Basic Language. The PPROBE program 
is shown in Appendix A. 

Manual control was used to initialize the pressure probe position 
prior to collecting data for the actual run, which utilized computer-controlled 
movements. While in the manual mode, the program asked a series of ques- 
tions which enabled the operator to test motor movements and set traverser 
motor velocity and acceleration default values. Through manual inputs, the 
probe was positioned such that the Pı (total pressure) hole was centered on 
the lengthwise axis of the VLSAM model body and placed as close as possible 
to it. From this point, the probe was moved vertically downward to the posi- 
tion of the first point in the data collection plane. For this experiment, the 
origin was located 2.75 inches below the model axis, and the data field dimen- 


sions were 2.75 inches above and below and 3 inches outward from the model. 
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After this initial probe position was set, the computer-controlled 
motor movement option was selected. The field dimensions (x, y coordinates), 
the traverser motor step distance, and the input file name were entered into 
the program. The data plane was 3.0" by 5.5" with a 0.25-inch step distance 
for this study. From this input, PPROBE then reiterated the total number of 
points to be measured (299 this case) and assigned filenames for each column 
of data. An example of how data files were named is as follows: 

example: ROA1A3 


where: 


run 
grid number (type) 

VLSAM model configuration 
field dimension (3, 5.5) 

step distance (0.25) 

test number 


OPr Pon 
ll H H H H H 


Thus, from the example above, the program assigned filenames 
R0A1A301.DAT through R0A1A313.DAT, which represented the 13 columns 
of data (23 points per column). 

The 5-hole pressure probe scale wheel was then adjusted until the 
Po and P3 (lateral) pressures were nearly equal (nulling), as measured by a 
portable digital manometer/calibrator. The measured yaw angle was read off 
the wheel and typed into the computer. Once this was accomplished, 
PPROBE moved the scanivalve from port 1 to port 4 via the Relay Actuator. 
There was a one-second delay to allow pressure equalization before the 
Digital Multimeter sampled the output voltage from the scanivalve trans- 
ducer via the Relay Multiplexer. After ten samples were taken at port 4, the 
Relay Actuator stepped the scanivalve to the next port (5), where another ten 
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samples were taken. This process was repeated until all five channel pres- 
sures (ports 4 through 8) were measured. [Ref. 3] Note that scanivalve port 4 
represents probe pressure P1, port 5 is P2, port 6 is P3, port 7 is P4, and port 8 

is P5. 

The Relay Actuator homed the scanivalve to port 48 after all the 
pressures were measured, and then PPROBE displayed the measurements 
and average values for each channel on the computer screen. The program 
either moved the probe upward one step (0.25") or remeasured the same 
point, depending on whether the data was within tolerance or not, as deter- 
mined by the operator. For this study, the following tolerances were used: P9 
and P3 differed by 0.09 psf or less, and P, was a positive number (or a very 
small negative number on the order of -0.5 psf). Once a column of data was 
measured, PPROBE would store the average values for each point in a file (23 
pts) and move the traverser to the next column position. The data acquisition 
process was continued until all 299 points were completed. 

2. CALP Program (Scanivalve Transducer Calibration) 

The CALP program (Appendix B) was the other data acquisition 
program utilized in this study. It was run both before and after the actual 
test (PPROBE) to account for any change in experimental conditions which 
might have occurred over the 8—10 hour period it took to run PPROBE. 

The transducer voltage was first adjusted to approximately zero 
millivolts. Calibration manometer (Figure 20) readings were then entered 
into the computer. The manometer provided a known pressure source for 


scanivalve calibration. From the transducer output voltage and pressure data 
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Figure 20. Calibration Manometer [Ref. 3] 


provided by CALP, two calibration curve equations were calculated (before/ 
after PPROBE runs). An averaged slope equation was then used in data 
reduction computations to minimize errors. This slope equation was used by 


the CONVERT program to change voltage data into dynamic pressure data. 
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3. CONVERT Program 
This program (Appendix C) was used to reduce data. It read the 
PPROBE data files (ex: ROA1A303.DAT) containing average pressures (Pl— 
P5) and converted them into x-y coordinates, velocity; yaw angle, pitch angle 
(alpha), total local pressure, total pressure coefficient, local static pressure, 
and static pressure coefficient. The x, y coordinates and yaw angle were sim- 
ply read from the PPROBE data files and input into CONVERT. The pitch 
angle coefficient was determined by a ratio of the pressures measured by the 
5-hole probe. This coefficient was then used in an equation developed by Lung 
[Ref. 3] using commercial curve-fitting software in order to find the corre- 
sponding pitch angle. Calibration curves provided by the probe manufacturer 
were used for the curve-fit. The program also obtained velocity pressure coef- 
ficients for particular pitch angles, from which the local velocity was calcu- 
lated after addition of a wind tunnel calibration factor, K, and other unit con- 
version factors. Additionally, the total and static pressure coefficients were 
obtained for different pitch angle regions. 
Room ambient pressure was used as the reference pressure for these 
coefficients, which are functions of total and static pressures and are non- 


dimensionalized by the tunnel dynamic pressure: 


Cps = (Psp — Ps)/Q (6) 
Cpr » (Bip - P4/Q (7) 
where 

Cps = Static pressure coefficient 

Cpr = Total pressure coefficient 

Q = Freestream dynamic pressure 

Ps = Freestream static pressure 

P. = Freestream total pressure 

Psp, = Local static pressure 

Pit = Local total pressure 


The actual dynamic pressure is (nearly) the same for the two cases. 
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For this experiment, the reference dynamic pressure values were 7.2 cm 
H20 for grid 0 and 10.0 cm H90 for grid 3. Temperature input for the 
CONVERT program was an average of the initial and final wind tunnel 
temperatures for the entire run time. Similarly, barometric pressure values 
were recorded before and after each run. CONVERT also added yaw and 
pitch angle (alpha) corrections to the data. These factors are +5.0° for yaw 
and —17.942° for pitch. They were determined from a preliminary run 
conducted with no grid and no missile model in the wind tunnel. Further 
explanation of this run is discussed in the Preliminary Tests section. The 
output of the CONVERT program was stored in a file named RESULT.DAT, 
which, in turn, was used as input to the TECPLOT system. 

4. TECPLOT 

The commercial TECPLOT software system was used to generate 
crossplane velocity vector plots and pressure contour plots. These plots could 
be tailored in many different ways by choosing scale factors, arrowhead 
wedge angles, contour levels and spacing, and many other parameters. A 
Hewlett-Packard 7470A x-y pen plotter was utilized in conjunction with 
TECPLOT to provide both the vector and contour plots. 


E. PRELIMINARY TESTS 
1. Dynamic Pressure Calibration 
All of the turbulence grids were previously calibrated in the wind 
tunnel. Readings from the tunnel calibration manometer and from a pitot- 
static tube inserted in the center of the test section were recorded over a 
speed range and wind tunnel calibration factors were obtained. These factors 


were used to adjust the tunnel flow velocity to the expected experimental 
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condition for the different grids. [Ref. 3] The calibration factors K are 0.8891, 
1.5084, 1.6487, 1.6545 and 1.1167 for no grid, and grids 1 through 4 respec- 
tively. These values were used in the CONVERT program to calculate the 
pressure and velocity in the test section. 

2. Yaw and Pitch Angle Corrections 

A test was conducted to find correction factors for yaw and pitch 
angle. From previous arrow plot data by Lung [Ref. 3], inconsistent crossflow 
velocity magnitudes and directions were noted toward the outer boundaries of 
the body-only missile configuration run. These outer boundaries represented 
the wind tunnel freestream region, where crossflow velocity is expected to 
reach zero. Therefore, to duplicate just the freestream region, this prelimi- 
nary test consisted of placing the pressure probe in the tunnel with no 
VLSAM model and no grid. Thus, the expected pitch and yaw angles should 
both have been zero. This was not the case however. 

Results of the preliminary run (R001A2) are listed in Appendix D 
(Result 00.DAT File). The yaw angles measured ranged from -4.00° to -7.00°, 
with an average of -5.00°. The pitch angles (alpha) ranged from approxi- 
mately +17.4° to approximately +18.4°, with the average +17.942°. Thus, to 
correct for these errors, +5.00° was added to the yaw angle and -17.942° was 
subtracted from the alpha values in the output file of the CONVERT program 
(RESULT.DAT). 

Though the exact cause for the yaw and pitch angle errors was not 
known, one possibility might have been a slight bend which was noted in the 
5-hole pressure probe. Other causes may have been improper calibration of 
the probe or a misalignment of the traverser assembly. The corrected errors 
only effect the velocity vector plots, and should have no effect on the pressure 


contour plots. 


4 


III. RESULTS 


The following sections discuss the velocity vector plots and the total and 
static pressure coefficient contour plots for the VLSAM model configurations 
A (plus) and C (cross), both with and without turbulence. All plots depict the 
3" by 5.5" data acquisition field and the position of the missile model (nose 
aspect) relative to the field. Vortex sizes and vortex distances from the model 


surface are referenced to the model base diameter d (1.75"). 


A. CONFIGURATION 0A ("PLUS" WITHOUT TURBULENCE) 

For the plus configuration, the swirling patterns of the velocity vector 
plot (Figure 21) clearly denote the two asymmetric vortices, which form 
circles that rotate in opposite directions. Although the bottom vortex center is 
evident, this plot fails to show the center of the top vortex. The vortex 
strength is a maximum on the outer edges of the vortex cores, denoted by the 
large vector arrows, where the velocities flow back toward the missile body. 
Towards the outer boundaries of the data acquisition field, where the vortex 
strengths are minimal, the vectors plot as points. 

The total pressure coefficient (CpT) contour plot (Figure 22) shows that 
the extent of the top vortex is approximately 0.72d at a distance of 0.33d from 
the missile surface. The bottom vortex is 0.83d at 0.67d from the body. There 


are more changes in the pressure gradient within a smaller area for the top 
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Figure 21. Velocity Vector Plot — Configuration 0A 
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Figure 22. Total Pressure Coefficient — Configuration 0A 
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vortex, which varies in Cpr from -0.14 to -1.35 within an inch, indicating 
greater vortex strength. The bottom vortex is slightly more diffused. The Cpr 
plot shows that the distance between the two vortex centers is approximately 
0.77d. 

The static pressure coefficient (Cps) contour plot (Figure 23) shows that 
the top vortex extends approximately 1.0d at a distance of 0.49d from the 
body. The bottom vortex extends about the same (1.1d) at a distance of 0.55d. 
Again, the top vortex appears to be stronger than the bottom vortex. Cps for 
the top vortex varies from -0.55 to -2.70 within 1.2 inches while the bottom 
vortex varies from -0.55 to about -2.18 in 1.6 inches. The distance between the 


vortices is approximately 0.70d on the Cps plot. 


B. CONFIGURATION 3A ("PLUS" WITH TURBULENCE) 

With added turbulence (grid 3), the velocity vector plot (Figure 24) still 
indicates vortex asymmetry, but it also indicates that the vortices have less 
strength (smaller vector arrows) than for the no grid condition. Again the 
vortex strength is maximized on the outer edges of the vortex cores. 

From the Cpr contour plot (Figure 25), the top vortex extends to approx- 
imately 0.72d, centered at about 0.44d from the model. The bottom vortex is 
slightly larger (0.83d) and is located 0.55d from the model body. The pressure 
gradient of the top vortex is steeper than for the bottom vortex, indicating 
greater strength. Cpr varies from -0.14 to -1.5 within an inch for the top 
vortex and from -0.14 to -1.35 within 1.4 inches for the bottom vortex. The 
distance between the two vortices as measured on the Cpr plot is approxi- 


mately 0.72d. 


= 
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Figure 23. Static Pressure Coefficient — Configuration 0A 
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Figure 24. Velocity Vector Plot — Configuration 3A 
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Figure 25. Total Pressure Coefficient — Configuration 3A 
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Examination of the Cps contour plot (Figure 26) reveals that both 
vortices are roughly 0.83d wide. The top vortex is closer to the missile body 
(0.39d) than the bottom vortex (0.55d). More contour levels in a smaller area 
and, consequently, a higher gradient and stronger vortex exists for the top 
vortex. Cps for the top vortex varies from -0.27 to -1.60 within 0.8 inch while 
the bottom vortex Cps ranges from -0.27 to -1.36 within 1.0 inch. On the Cps 


plot, the two vortices are about 0.66d apart. 


C. CONFIGURATION 0C ('CROSS" WITHOUT TURBULENCE) 

For the model C configuration (cross), the velocity vector plot (Figure 27) 
still clearly displays the two asymmetric vortices. As with configuration A, 
the strengths of the vortices are maximized at the edges of the circular swirls, 
where the flow is back toward the missile body. 

Figure 28, the total pressure coefficient plot, displays a 0.77d wide top 
vortex located at about 0.5d from the model surface. The bottom vortex 
extends to approximately 0.88d and is 0.66d from the model. The bottom 
vortex is more tightly wrapped (i.e., more contour levels per area) closer to its 
core than is the top vortex. The Cpr range for the bottom vortex (-0.14 to 
-1.23) 1s slightly greater than the range for the top vortex (-0.14 to -1.09). The 
vortex centers are roughly 0.66d apart. 

From the static pressure coefficient contour (Figure 29), both vortices 
extend to 0.95d with the top vortex slightly closer to the model body (0.56d) 
than the bottom vortex (0.66d). The relative strengths of the vortices is diffi- 
cult to interpret from the Cps plot, but the top vortex appears to be a little 
less diffused than the bottom vortex and therefore stronger. As with the total 


pressure coefficient plot, t::2 vortices are 0.66d apart. 


92 


Y-AXIS (IN) 





1.0 ου 


0.0 


-2.0 -1.0 
X-AKIS (IN) 


Figure 26. Static Pressure Coefficient — Configuration 3A 
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Figure 27. Velocity Vector Plot — Configuration 0C 
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Figure 28. Total Pressure Coefficient — Configuration 0C 
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Figure 29. Static Pressure Coefficient - Configuration 0C 
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D. CONFIGURATION 3C (‘CROSS" WITH TURBULENCE) 

The velocity vector plot (Figure 30) shows that the addition of turbulence 
(grid 3) diffuses the asymmetric vortices for model configuration C. The vector 
arrows plot smaller for this turbulent case than for the no grid case. 

The Cpr contour plot (Figure 31) reveals a 0.92d bottom vortex at a dis- 
tance of 0.66d from the model body and a 0.83d wide top vortex at 0.55d from 
the model. The bottom vortex appears to be slightly more tightly wrapped 
towards the vortex center than the top vortex. Cpr levels for both vortices 
range from about 0.0 to -1.23. The distance between the vortices is approxi- 
mately 0.72d. 

Figure 32, the static pressure coefficient plot, like the Cpr plot, also 
shows the distance between the vortices to be 0.72d. The top vortex (extent 
0.77d) is roughly 0.55d from the model's surface while the bottom vortex 
(extend 0.88d) is 0.66d. The vortices appear to have the same strength on the 


Cps plot, in which the coefficients range from 0.0 to approximately -1.1. 


E. COMPARISONS 
1. Between Body Configurations (A and C) 
The following observations were noted when making comparisons 
between the two missile configurations for both the turbulent (grid 3) and 
non-turbulent (no grid) runs. 


(1) From the total pressure coefficient (Cpr) contour plots, the vortices 
are slightly larger for configuration C. This is more pronounced for 
the run conducted with added turbulence (grid 3). 
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Figure 30. Velocity Vector Plot — Configuration 3C 
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Figure 31. Total Pressure Coefficient — Configuration 3C 
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Figure 32. Static Pressure Coefficient — Configuration 3C 
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(2) 


(3) 


(4) 


2. 


The vortices of the C configuration appear to be centered a little fur- 
ther away from the model surface than do the A configuration 
vortices. This shift is on the order of 0.2 inches, and is apparent in 
both the no-grid and grid 3 runs. 


The relative distance between vortex centers does not seem to be 
affected when the configurations are changed. 


The bottom vortex is more tightly wrapped (stronger) for configura- 
tion C based on the Cpr plots. This is more noticeable for the run 
conducted without turbulence. 


Between Turbulence Levels (0 and 3) 


The following observations were noted when making comparisons 


between the two turbulence levels tested for each body configuration. 


Turbulence level 0 was the no-grid condition and level 3 was the run con- 


ducted with a grid inserted. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


The velocity vector plots indicate smaller arrows, therefore weaker 
vortices (less crossflow), when turbulence is added. This trend holds 
for both body configurations. 


The two vortices remain in approximately the same positions as indi- 
cated on both the vector plots and the static pressure coefficient (Cps) 
plots). 


From the total press coefficient (Cpr) plots, the center of the bottom 
vortex for configuration A shifts slightly closer to the model body 
when the turbulence increases. (-0.2 to 0.0 on the x—axis). 


Also on the Cpr plots, with added turbulence, the top vortex core for 
configuration C seems to shift away from the VLSAM body. 


The Cpr plots for both configurations are slightly more diffused with 
added turbulence. Specifically, the vortex centers are less tightly 
wrapped, thus indicating weaker vortices for the more turbulent 
condition. 


There are no noticeable differences in the static pressure coefficient 
contours between the turbulent and non-turbulent runs. 
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3. With Body-Only Configuration (B) 

The following observations were made in order to correlate the 
results of this study with the previous experiment conducted by Lung for a 
body-only VLSAM model configuration. [Ref. 3] The purpose is to note how 
the addition of wings and strakes might affect the vortex flow pattern around 
the missile. The velocity vector, total pressure and static pressure coefficient 
plots for the body-only run (0B) are displayed in Figures 33, 34, and 35 
respectively. These comparisons were made for the nominal ambient (no-grid) 
flowfield condition. 


(1) From the velocity vector plots, the vortex pattern for configuration B 
is similar to the one for body C in that the vortex asymmetry is more 
pronounced for these cases. Body A vortices are closer to the model 
surface and asymmetry is less pronounced. 


(2) From the total pressure coefficient contours, the vortices appear to be 
located much closer to the missile body for body A than for either 
body C or B. Apparently the strake/wing-generated vortices for this 
configuration act to hold the nose-generated vortices near the body, 
resulting in the higher induced side forces for this configuration 
observed by Rabang. 


(3) Also from the CpT plots, the strengths of the bottom vortices appear 
to be roughly the same for bodies A and B, while the plot for body C 
seems to be more tightly wrapped (stronger) near its center. 


(4) From the static pressure contours, the body-only vortices are more 
diffused and larger than either configuration with wings (A and C). 


(5) The relative distances between the top and bottom vortices remains 
the same on all three plots for all body configurations. 


In general, the vortex pattern around the nose of the missile model with 
wings (in either configuration) resembles the vortex pattern for a body-only 


configuration. Though there are subtle differences as previously noted, the 
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relative strengths, sizes and positions of the asymmetric vortices were com- 


parable in all cases. 
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Figure 33. Velocity Vector Plot — Configuration 0B 
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Figure 34. Total Pressure Coefficient — Configuration 0B 
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Figure 35. Static Pressure Coefficient — Configuration 0B 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The flowfield about a Vertically-Launched Surface-to-Air Missile 


(VLSAM) model at high angle of attack was investigated in the wind tunnel 


of the Naval Postgraduate School. Missile "plus" and "cross" body configura- 


tions (A and C respectively) were both tested. The angle of attack was set at 
50° with a Reynolds number of 1.1x10° for all data runs.Two flowfield condi- 


tions were treated: the nominal ambient wind tunnel condition (no grid) and 


a condition with a grid-generated turbulence of length scale 1.08 inches and 


1.88% turbulence intensity (see Table 2). The following conclusions were 


reached: 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


An increase in turbulence intensity tended to reduce the strength of 
the asymmetric nose-generated vortices. 


The two asymmetric vortices remained in approximately the same 
position for an increase in turbulence. 


The top vortex was closer to the model surface and appeared to be 
stronger for both body configurations. This condition was more pro- 
nounced for configuration A. 


The wing/strake arrangement of configuration C caused the vortices 
to be centered further away from the model surface than those of con- 
figuration A, correlating with the differences in induced side forces 
for those configurations observed by Rabang. 


Configuration C vortices are more diffused and larger. This was more 
apparent when turbulence was added. 


Crossflow velocity vector plots agreed with the behavior denoted by 


the total and static pressure coefficient contours for both body config- 
urations and turbulence levels. 
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(7) Though subtle differences exist, the addition of wings and tails did 
not greatly alter the vortex pattern around the nose of the missile 
model, when compared to the body-only configuration tested by Lung. 


Recommendations for a continued study of the behavior of asymmetric 

vortices under varying flowfield conditions are suggested as follows: 
(1) Examine the vortices at positions further back along the model body, 
where effects from the wings might better be seen. (Such as at a 


length/diameter ratio of 9.) 


(2) Investigate asymmetric vortex behavior for just a body-only configu- 
ration, at the position described above, to provide comparisons. 


(3) Continue to study vortex behavior at various angles of attack and 
turbulence with varying intensities and length scales in order to pro- 
vide a large data base, which can be used to calculate vorticity 
contours. 
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APPENDIX A. PPROBE PROGRAM 


DEF SEG:CLEAR ,&HFEO0:GOTO L 'Begin PCIB Program Shell 
GOTO 1000 ' User program 

GOTO 900 ' Error handling 

I-&HFEOO ' Copyright Hewlett-Packard 1984,1985 
PCIB.DIR$-ENVIRON$("PCIB") 

I$T-PCIB.DIR$* 'NPCIBILC.BLD" 

BLOAD I$,I 

CALL I(PCIB.DIR$,1%,J%):PCIB.SEG=1% 

IF J%=0 THEN GOTO 13 

PRINT "Unable to load."; 

PRINT " (Error hija, ΠΠ; 

END 


DEF SEG=PCIB.SES:O.S=5:G.S=106 1° U το 
J.G 2ooO:D p 25:LDD EIDE 30 
GET .MEM=35:L.S=40:PANELS=45:DEF.ERR=50 


PCIB.ERR$-STRING$(64,32) : PCIB.NAMES=STRING$(16,32) 

CALL DEF.ERR(PCIB.ERR,PCIB.ERR$,PCIB.NAMES,PCIB.GLBERR) : PCIB.BASERR-255 
ON ERROR GOTO 3 

J=-1 


I$=PCJB.DIR$+"\PCIB. SYN" 

CALL O.SCIS$) 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

1:0) . 

CALL I.V(1,READ.REGISTER, READ. SELFID, DEFINE, INITIALIZE ου το αν, 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL 1.V(1,ENABLE.SYSTEM, DISARLE. SYSTEM, INITIALIZE, POWER ος, 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL I.V(I, MEASURE OUTPUF START HAEE) i 
IF PCIB.ERR<>0 THEN ERROR PCIB. BASEFRR 

CALL 1.V(1,ENABLE. INT.TRIGGER, DISABLE. INT. TRIGGER, ENABLE. OUTPUT, DISABLE. OUTPU 


IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 

CALL I.V(I,CHECK.DONE,GET.STATUS SET-BUBNGTTIONOSEROPANSES 
IF PCIB.ERR<>0 THEM ERROR, PCIB. BASERR 

CALL 1.V(1,SEP MODE, WRITE SCA READ CAL, cTORE. CAL 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 

CALL 1.V(1,DELAY SAVE. SYSTEM ΠΠ 

IF PCIB-ERR<>O0 THEN ERROR PCIB.BASERR 

I-1 

CALL I.V(1I1,SET.GATETIME,SET.SAMPLES, SET SLOPE, SET SOUERECEJ 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

Շուն I.C{(1,FREOUENCY , AUTO. FREQ, PERIOD AUTOTPER) 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL I.C(I,INTERVAL,RATIO,TOTALIZE,R100MILLI) 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL I.C(I,R1,R10,R100,R1KILO) 

IF PCIB.ERR<K>O THEN ERROR PCIB.BASERR 

CALL I.C(I,R10MEGA,R100MEGA, CHAN.A,CHAN.B) 

IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

CALL I.C(I,POSITIVE,NEGATIVE,COMN, SEPARATE) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 

152 

1-3 
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CALL I.V(I,ZERO.OHMS, SET.SPEED,J,J) 
IF PCIB.ERR<>0 ‘THEN ERROR PCIB.BASERR 
CALL I.C(I,DCVOLTS,ACVOLTS,OHMS,R200MILLI) 
IF PCIB.ERR«&?»0 THEN ERROR PCIB.BASERR 
CALI. I.C( I,R2,R20,R200,R2KILO) 
IF PCIB.ERR<>O THEN ERROR FCIB.BASERR 
CALL I.C(I,R20KILO,R200KILO, R2MEGA,R20MEGA) 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
CALL 1. CCl; AUIOM, R2.5,R1275,.)) 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
I=4 
CALL I.V(I,SET.COMPLEMENT,SET.DRIVER,OUTPUT.NO.WAIT,ENABLE.HANDSHAKE) 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
CALL I.V(I,DISABLE.HANDSHAKE,SET.THRESHOLD,SET.START.BIT,SET.NUM.BITS) 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
CALE I. V(I, SET.LOGIC.SENSE, J, J.J) 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
CALL I.C(I,POSITIVE,NEGATIVE, TWOS,UNSIGNED) 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
CALE). CUl Oc, LIL. RO.RI) | 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
CALL I.C(I.R2,R3.R4, R5) 
ΙΕ PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
ΙΙΙ Γ᾽ ΠΠ5) Հ 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
ΓΗ ΙΓ RIO RII. RI2;R13} 
IE FCIB.ERR<>20 THEN ERROR. PCIB.BASERR 
CALL I.C(I,RIS,R15.R16,J) 
IF PCIB.ERRK20 THEN ERROR FPCIB.BASERR 
I*6 
Ge VUISSET.FREQUENCY SET C AMPLITUDESSET.OPEFSETI.SET.SYMMETRY) 
IE PCIB.ERR<>0 THEN ERROR PCIB. BASERR i 
Տուն Vil Sel. BURST COUNT J ,J , T) 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
ΗΠ ΤΕΙ clLNE, SOUARP  JRIANGEE, CONTINUOUS } 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
ΓΗ SO (I GATED: BURSE, Jd) 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
177 
CIHEBSDOVODSADUIOSCRHERSOCADIBRADES SETSSEHMSITIVITY SET.VERITI.OFFSET) 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
CALCI VI, SET. COUFLING, SET POLARITY: SET.: SWEEFSPEED,SET.DELRY)? 
IF PCIB.ERR<>O THEN ERROR PC] B.BASERR 
CALLVI-V (1, SET. TRIG SOURCE SET.TRIG:SEOPE,SET.TRIG.LEVEL,SET.TRIG.MODE) 
IF PCIB.ERR<>?0 THEN ERROR PCIE. BASERR 
CALL Ἱ νι ους ΟΘΙΝΟΕΕ WF GET.TWO.WE GET.VERT.INEFO,GET.TIMEPASE. IHFO) 
IF PCIB.ERR<>0 THEN ERROR FCIB.BASERR 
CALL I. V(I, GET.TRIG.INFO,CALC.WFVOLT,CALC.WFTIME,CALC.WF.STATS) 
TE PCIR: ERR<>0O THEN ERROR PCIB.BASERR 
CALL I. V(I. .CALC.RISETIME,CALC.FALLTIME,CALC.PERIOD,;CALC. FREQUENCY ) 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
CALL 1.V(1,CALC.PLUSWIDTH, CALC.MINUSWIDTH, CALC.OVERSHOOT , CALC. PRESHOOT ) 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
CREE we Vr, CALC. ER. TO ΕΡΕ ΞΕ; Πο SCOPE -SINRT, MEASURE. SINGLE. WE) 
IF PCIB.ERR<>0 THEN ERROR PCIB. BASERR 
CALL I.V(I,MEASURE.TWO.WF,J,J,J) 
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IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
CALL J.C(I,RIONANO,R100NANO,R1MICHRO,R10MICRO) 
IF PCIB.ERR<>0 TIIEN FRROR PCIB.BASERR 


CALL 1.C(1,R1IOOMICRO,RIMILLI,R1OMILLI,RIOOMILLI ) 


IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
CALL I.C(I,R1,R30,R20NANO,R200NANO) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 
CALL I.C(I,R2MICRO,R20MICRO,R200MICRO,R2MILLI) 
IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 
CALL I.C(I,R20MILLI,R200MILLI,R2,R20) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 
CALL I.C(I,R50NANO,R500NANO,R5MICRO,R50MICRO) 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
CALL 1.C(1,RSO0OMICRO,RSMILLI,RSOMILLI,RSOOMILLI ) 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
CALL 1I.CC1,R5, R50, CHAN GA. CHAN Է) 

IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
CALL 1I.C(1, EXTERNAL, POSITIVE, NEGATIVE TAC) 
IF PCIB.ERR<>O THEN ERROR PCIB.RËSERR 
CALL I.C(I,DC,TRI GGERED; AUTO. TRIG AUTOSBEVED) 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
CALL I.C(I,X1,X10,STANDARD, AVERAGE) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 

1-8 

CALL 1.V(1],Q0PEN. CHANNEL, CLOSE -CHANNEI. 0.0? 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
CALE CTS š 

IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
IS=PCIB.DIR&S+'"\PCIB. PID" 

CALL L.P(I$) 

IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
If="DMM.01"3:1=32:J=0:K=0:L=1 

CALL DEFINE COMM. Oils ln kL 

IF PCIB.ERR<>O THEN ERROR PCIR.BASERR 
I$="Func.Gen.01":1=6 J 0 Kk j P] i 
CALL DEFINE(Func.Gen.01,I$,I,J,K,L) 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
Itft="Scope.01":I=7:J=0:K=2:L=1 

CALL DEFINE(Scope.01,1$,1,J,K,L) 

IF PCIB.ERR<>O THEN ERROR PCJB.BASERR 
yts"Counter.01":1-1:J-*0:K*3:L*1 

CALL DEFINE (@ounter 01, 12.1 JOE) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 
If="Dig.In.01":1=4:J=0:K=4:L=1 

CALL DEFINE(CDYg.In.0DI P$ PO 1) 

IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
If="Dig.Out.01":1=4:J=1:K=4:L=1 

CALL DEFINE( Die [Out. 01, 14.1.2...) 

IF PCIR.ERR<>O THEN ERROR PCIB.BASERR 
I1$="Relay .Act.01":1=8:J=0:K=5:L=1 

CALL DEFINE(Relay.Act.01,1$,1,J,K,L) 

IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
I$="Relay .Mux.01":1=2=2:J=0:K=6:L=1 

CALL DEFINE(Relay.Mux.01,1$,1,J,K,L) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 
IS=ENVIRON$ ( "PANELS" )+'"\PANELS.EXE" 

CALL L.S(1I$) 

GOTO 2 
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900 IF ERR=PCIB.BASERR THEN GOTO 903 

901 PRINT "BASIC error #":ERR:" occurred in line ":ERL 

902 STOP 

903 TMPERR=PCIB.ERR: IF TMPERR=0 THEN TMFERR=PCIB.GI.BERR 

906 PRINT "PC Instrument error W&";TMPERR;" detected at line ";ERL 

905 PRINT “Error: ";PCIB.ERR$ 

906 IF LEFT$(PCIB.NAME$,1)<>CHR$(32) THEN PRINT “Instrument: ";PCIB.NAMES 

907 STOP 

908 COMMON PCIB.DIR$,PCIB.SEG 

909 COMMON LD.FILE,GET.MEM, PANELS, DEF.ERR 

910 COMMON PCIB.BASERR,PCIB.ERR, PCIB.ERR$€, PCIB.NAME$, PCIB.GLBERR 

911 COMMON READ.REGISTER,READ.SELFID,DEFINE, INITIALIZE.SYSTEM, ENABLE.SYSTEM, DISA 
BLE .SYSTEM, INITIALIZE, POWER.ON, MEASURE, OUTPUT ,START,HALT, ENABLE. INT.TRIGGER,DISA 
BLE. INT.TRIGGER, ENABLE.OUTPUT, DISABLE.OUTPUT, CHECK. DONE, GET. STATUS 

912 COMMON SET.FUNCTION,SET.RANGE,SET.MODE,WRITE.CAL,READ.CAL,STORE.CAL,DELAY,SA 
VE.SYSTEM,SET.GATETIME,SET.SAMPLES,SET.SLOPE, SET. SOURCE, ZERO.OHMS,SET.SPEED,SET. 
COMPLEMENT,SET.DRIVER,OUTPUT.NO.WAIT,ENABLE.HANDSHAKE,DISABLE.HANDSHAKE 

913 COMMON SET.THRESHOLD,SET.START.BIT,SET.NUM.BITS,SET.LOGIC.SENSE, SET. FREQUENC 
rose, AMPLITUDE, SET OFFSET ,SEIT.SYMMETRY,SET. BURST. COUNT, AUTOSCALE,CALIBRATE,SET. 
BENSTIIVITY,SET.VERT.OFFSET,SET.COUPLING,SET.POLARITY,SET.SWEEPSPEED 

914 COMMON SET.DELAY,SET.TRIG.SOURCE,SET.TRIG.SLOPE,SET.TRIG.LEVEL,SET.TRIG.MODE 
,GET.SINGLE.WF,GET.TWO.WF,GET.VERT. INFO, GET. TIMEBASE.INFO,GET.TRIG.INFO,CALC.WFV 
OET,CALC.WFTIME,CALC.WF.STATS,CALC.RISETIME,CALC.FALLTIME,CALC.PERIOD 

915 COMMON CALC.FREQUENCY,CALC.PLUSWIDTH,CALC.MINUSWIDTH, CALC. OVERSHOOT, CALC.PRE 
SHOOT,CALC.PK.TO.PK,SET.TIMEOUT,SCOPE.START,MEASURE.SINGLE.WF,MEASURE.TWO.WF,OPE 
N.CHANNEL, CLOSE. CHANNEL 

916 COMMON FREQUENCY,AUTO.FREQ,FERIOD,AUTO.PER,INTERVAL,RATIO, TOTALIZE,R100MILLI 
,R1,R10,R100,R1KILO,R10MEGA,R100MEGA,CHAN.A,CHAN.B, POSITIVE,NEGATIVE, COMH, SEPARA 
TESDCVOLTS,ACVOLTS,OHMS,R200MILLI,R2,R20,R200,R2KILO,R20KILO,R200KILO 

917 COMMON R2MEGA,R20MEGA,AUTOM,R2.5,R12.5,POSITIVE, NEGATIVE, TWOS,UHSIGNED,OC, TT 
L,RO,R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11,R12,R13,R15,R15,R16,SIHE, SQUARE, TRIANGLE 
,CONTINUOUS, GATED, BURST,R10NANO,R100NANO,R1MICRO,R10MICRO,R100MICRO 

918 COMMON RI1MILLI,R10MILLI,R100MILLI,R1,R10,R20NAHO,R200NAHO,R2MICRO,R20MICRO,R 
200MICRO,R2MILLI,R20MILLI,R200MILLI,R2,R20,R50HAIHO,R500HANO, RSMICRO,R50MICRO,R50 
OMICRO,RSMILLI,RSOMILLI,RSOOMILLI,R5,R50,CHAN.A,CHAN.B, EXTERNAL, POSITIVE 

919 COMMON NEGATIVE,AC,DC, TRIGGERED, AUTO. TRIG, AUTO. LEVEL, X1,X10, STANDARD, AVERAGE 
920 COMMON DMM.01,Func.Gen.01,Scope.01,Counter.01,Dig.In.01,Dig.Out.Ol.Relay.Act 
.01,Relay.Mux.01 

999 'End PCIB Program Shell 

1000 REM This step initialzes the HP system 

Lon CLS 

ο ο ΟΕΤΙΟΝ BASE 1 

մլ ինն ի P(5),PA(50,5),.PP(50,5),XPT(50),YPT(50) ,X(50),Y(50) , YAW(50) 

1040 REM 

1050 CALL INITIALIZE.SYSTEM( PGMSHEL.HPC) 

1060 REM 

1070 REM SET FUNCTIOHOH THE 'DMM' , 'RELAY MUX , 'RELAY ACTUATOR' 

1080 REM 

1090 CALL SET. FUNCTION(DHM.01, DCVOLTS ) 

1100 CALL SET.RANGE(DMM.01, AUTOM) 

1110 CALL DISABLE.INT.TRIGGER(DMM.01) 

1120 CALL ENABLE.OUTPUT(RELAY.MUX.01) 

1130 CALL ENABLE.OUTPUT(RELAY.ACT.01) 

1140 REM ************* PROGRAM TRAVERSE *****x****3x**x*kx*xtx 

1150 REM 

1160 REM OPEN THE COM PORT AND INITIALIZE THE MOTOR SETTINGS 

πα. ο ΟΡΕΝ "comi:1200,n,8,1,rs,cs,ds,cd" AS #1 

1180 REM SET MOTOR DEFAULT VALUES 

MEIO DATA 2000,2000,2000,2,2,2,0.000125,0.000125,0.000125 , 
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READ V1,V2,V3 ,R1 R2 R3 


σα το 3 


REM DEFINE CHARACTERS FOR DATA REDUCTION ALGORITHM 
RN2$-"RENAME A:RAW.DAT " 


HEAD1$ - " 8 X x P1 P2 P3 P4 P5 YA 
FORMATS= "## 8.8 an. A8 WAH. HHRH AAR. HHR RN. ΝΗ. BHR ARA AAN 
PRINT 

PRINT "I OR ior ik kk kok 

PRINT "** USER MUST SELECT 'CAPS LOCK' FUNCTION **" 

PRINT PXXHKEKËKKEHEKËKEKKKKIEKKKEKËKKKETKKKKKKEKAN Y 

REM DISPLAY MOTOR DEFAULT SETTINGS 

PRINT " JOO Ok ROR eoe ob ooo oec oe o v e e oe n 

PRINT " INITIALIZED VALUES FOR ALL MOTOR SETTINGS:" 

PRINT " VELOCITY = 1000 STEPSZSEGI 

PRINT '' RAMP(MOTOR ACCELERATION) = 2 (6000 STEPS/SEC^2)" 
PRINT Ὁ DEFAULT INCREMENTAL UNITS ARE INCHES" 

PRINT " Xo oko ook koe ok oko doo X oe eX oe oe ooo boe oe OG o x 

PRINT 

PRINT "NOTE! ] USE MANUAL CONTROL TO INITIALIZE PROBE POSITION BEFORE" 
PRINT SELECTING COMPUTER CONTROLLED MOVEMENT. E 

PRINT 

INPUT "MANUAL CONTROL OR COMPUTER CONTROL (ENTER 'MAN' or 'CP')";CON$ 


IF CONS="CP" 


THEN 3490 


REM OPTION TO CHANGE DEFAULT SETTINGS OF VELOCITY OR ACCELERATION RAMP 


PRINT ' 
PRINT 

PRINT " DO YOU WANT TO CHANGE THE VELOCITY OR ACCELERATION RAMP" 

PRINT " DEFAULT SETTINGS? (Y or N)" 

PRINT : 

PRINT "IF 'NO', THIS PROGRAM WILL THEN LET YOU DEFINE THE" 

PRINT "DISTANCE YOU WANT TO MOVE (IN INCHES). IF 'YES'," 

PRINT "YOU CAN CHANGE ANY OR ALL OF THE DEFAULT SETTINGS FOR ANY MOTOR." 
PRINT 
PRINT i 
PRINT 

INPUT "DO YOU WANT TO CHANGE ANY OF THE DEFAULT SETTINGS? (Y or N)";D$ 
IF D$-"y" THEN 1590 

IF D$-"N" THEN 2220 

REM 

REM **** OPERATOR SELECTED MOTOR VARIABLES ***** 

PRINT 

PRINT 

INPUT “WHICH DEFAULT VALUE? (ENTER '1'FOR VELOC OR '2' FOR ACCEL RAMP)":L 
ON L GOTO 1690,1930 

PRINT "DO YOU WANT TO CHANGE THE DEFAULT VELOCITY? (Y OR N)" 

INPUT V$ 

IF V$-"y" THEN 1690 

PRINT "DO YOU WANT TO CHANGE THE DEFAULT ACCELERATION RAMP? (Y or N)" 
IF R$="Y" THEN 1990 

IF R$-"N" THEN 1450 

PRINT 

PRINT 

INPUT "WHICH MOTOR VELOCITY DO YOU WISH TO CHANGE? (1,2. or 3)";J 

ON J GOTO 1730,1830,1880 

PRINT 

PRINT 

INPUT "ENTER DESIRED VELOCITY OF MOTOR 51";V1 

PRINT 

PRINT 

PRINT 

PRINT "DO YOU WANT TO CHANGE VELOCITY OF ANOTHER MOTOR? 


(Y OR N)" 
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1800 INPUT v$ 

1810 IF V$-"Y" THEN 1690 

1820 IF V$="N" THEN 1430 

1830 PRINT 

1840 PRINT 

NSSOCINPUT "ENTER DESIRED VELOCITY OF MOTOR 2";V2 
1860 PRINT 

1870 GOTO 1780 

1880 PRINT 


1890 PRINT 

1900 INPUT "ENTER DESIRED VELOCITY OF MOTOR $3";V3 
1910 PRINT 

1920 GOTO 1780 

1930 PRINT 


1940 PRINT 

1950 INPUT "WHICH MOTOR ACCEL RAMP DO YOU WANT TO CHANGE? (1, 2, or 3)':K 
1960 ON K GOTO 1970,2060,2120 

1970 PRINT 

1980 PRINT 

1990 INPUT "ENTER DESIRED ACCELERATION RAMP OF MOTOR 41";R1 

2000 PRINT 

2010 PRINT 

2020 PRINT "DO YOU WANT TO CHANGE THE ACCEL RAMP OF ANOTHER MOTOR? (Y or N)?" 
2030 INPUT RM$ 

2040 IF RM$-"y" THEN 1930 

2050 IF RM$="N" THEN 1450 

2060 PRINT ` 

2070 PRINT i 
2080 INPUT “ENTER DESIRED ACCELERATION RAMP OF MOTOR 82':R2 

2090 PRINT 

2100 PRINT 

2110 GOTO 2000 

2120 PRINT 

2130 PRINT 

2140 INPUT "ENTER DESIRED ACCELERATION RAMP OF MOTOR 43";R3 
2150 PRINT 

2160 PRINT 

2170 GOTO 2000 i 

2180 REM 

2190 REM DEFINE DISTANCE TO MOVE MOTOR 

2200 PRINT 

2210 PRINT 


2220 PRINT 

2230 REM INITIALIZE MOTOR INCREMENTS TO ZERO 

2240 I1-0 

2250. 1250 

2260 1350 

2270 PRINT 

2280 PRINT " **»***»*****x**&t£t tx tttttttk*xttkitivkx*xxt*vxtxkxttttxktitixrxwrgn 
2220 PRINT " ** DEFINE WHICH MOTOR YOU WANT TO MOVE ας 
2300 PRINT " ** rx 
2310 PRINT " ** NOTEM IA POSITIVE (t) INCREMEUT TO-A MOTOR er 
2320 PRINT " ** MOVES TRAVERSER AWAY FROM THAT PARTICULAR MOTOR ru 
2330 PRINT " ** ee 
νυ PRINT * ee -- MOTOR fl MOVES THE PROBE UPSTREAM AGAINST THE FLOW **" 
2350 PRINT " ** -- MOTOR #2 MOVES THE PROBE TOWARD THE ACCESS WINDOW τ 
2360 PRINT " ** -- MOTOR {3 MOVES THE PROBE VERTICALLY DOWNWARD ΕΤ 
νκγ5.ννν------------Ξ---;--΄-.---“”--------΄-------“-----μπι 
2380 PRINT 

2390 PRINT 


2400 INPUT "WHICH MOTOR DO YOU WANT TO MOVE? (1,2, or 3)";L 
2410 ON L GOTO 2420,2680,2970 V 
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2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2910 
2920 
2930 
2540 
2990 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 
3070 


PRINT 

PRINT 

PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #1?" 

PRINT " ***x***** (ENTER DISTANCE IN @INCHES )” ** 22 eee 
INPUT I1 


PRINT 

PRINT" KEKAKKKKKKKKKKEKKKKKKKEKKKKKKKK KKK i 

PRINT 

PRINT "SUMMARY OF OPERATOR INPUTS:" 

PRINT " MOTOR 11 VELOCITY =V] 

PRINT " ACCELERATION RAMP = ";R1 

PRINT " INCREMENTAL DISTANCE = ";11;"INCHES" 


PRINT" *'* oe eb ne x x oe s ox o eo o eco x oko eoo KKK KEKE 


PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or N)" 
PRINT 

PRINT "ENTER 'N' TO START MOTOR MOVEMENT. ENTER 'Y' TO RETURN" 
PRINT "TO VARIABLE SELECTION SUBROUTINE." 

INPUT V$ . 

IF V$-"Y" THEN 1430 

GOSUB 3410 

PRINT 

PRINT "DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? (Y or N)?" 
INPUT C$ w 
ΙΕ ΟΥ THEN 2220 

IF C$="N" THEN 3260 

PRINT ° 

PRINT 

PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #2?" 

PRINT " ***se*t*** (ENTER DISTANCE IN ΠΕΠΕΞΡΣΣΣΙ ΤΑ 
INPUT I2 

PRINT 

PRINT 


REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 
PRINT" ΧΚΧΧΧΧΧΧΧΚΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ!Ι! 


PRINT 

PRINT "SUMMARY OF OPERATOR INPUTS:" 

PRINT " MOTOR 82 VELOCITY =) 3V2 

PRINT " ACCELERATION RAMP = ";R2' 

PRINT " INCREMENTAL DISTANCE = ";I2;"INCHES" 
PRINT" *kkkkkkkkkkkkkkkkkkkkkkkkkkkkkk*kx" 

PRINT 

PRINT 

PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or N)" 
PRINT 


PRINT "ENTER 'N' TO START MOTOR MOVEMENT. ENTER 'Y' TO RETURN" 
PRINT "TO VARIABLE SELECTION SUBROUTINE." 

INPUT V$ 

IF V$="Y" THEN 1430 

GOSUB 3410 

PRINT 

PRINT "DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? (Y or N)?" 
INPUT C$ 

IF C$="Y" THEN 2220 

IF C$="N" THEN 3260 

PRINT 

PRINT 

PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR $3?" 

PRINT " ********* (ENTER DISTANCE IN INCHES) ****xsxxx" 
INPUT I3 

PRINT 

PRINT 

REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 
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3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3339 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3690 
3660 
3670 
3680 
3690 
3700 


PRINT" ***-* ox x x xx xo ooo oO x x x eo X x e X PG x"! 


PRINT 

PRINT "SUMMARY OF OPERATOR INPUTS:" 

PRINT ” MOTOR 83 VELOCITY = ::03 

PRINT " ACCELERATION RAMP = ";R3 

PRIHT " SINCREMEHIAL DISTANCE - "T3 INCHES" 
PRINT 

PRINT" *xxoookokokoko oec ok x oc e o oc e o kc e e xo e X € € I e e" 

PRINT 

PRINT 

PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or N)" 


PRINT 

PRINT "ENTER 'N' TO START MOTOR MOVEMENT. ENTER 'Y' TO RETURN" 

PRINT "TO VARIABLE SELECTION SUBROUTINE." 

INPUT v$ 

IF V$="Y" THEN 1430 

GOSUB 3410 

PRINT | 
PRINT ! 
INPUT "DO YOU WANT TO INPUT ANOTHER MANUAL MOTOR MOVEMENT (Y or N)";M$ 
IF M$-"Y" THEN 2210 


PRINT 
PRINT "DO YOU WANT TO INPUT COMPUTER CONTROLLED MOTOR MOVEMENT?" 
PRINT " ********* NOTE!!! f****x*x*x*x*xx* "Ո" 


PRINT " ALL PREVIOUS MOTOR INCREMENT INPUTS HAVE BEEN ZEROIZED." 
PRINT "PROGAM WILL LET YOU CHOOSE MANUAL OR CP-CONTROLLED MOVEMENT." 
PRINT "***** (IF 'NO', THE PROGRAM WILL END). *****" 

PRINT 

INPUT "DO YOU WANT COMPUTER CONTROLLED MOTOR MOVEMENT (Y or N)";N$ 
IF N$2"Y" THEN 3500 


PRINT 

PRINT 

PRINT 

PRINT " SOOO OO c XO XX X X X X X 

PRINT " THE PROGRAM HAS ENDED." 

PRINT " Ֆո փի փին իի ԱԵՒ Ւ 

END 

REM ******* MOTOR MOVEMENT SUBROUTINE ΣΑΧΑΧΧΣΧΧΧ 

ԼԱ Ա oo PRANONI eMC Tas Cie oC 2a" 2 Ole": C3="3 C3 
PRINTS EP PI 1=Sus νι; RISI, ' 
ΕΠ |2=" το ο Om VA R27" R2 

PRINT 41, "I3s":I3;":V35s";V3:":R35s";R3;":Q" 

RETURN 

REM £ £ £ X Kk £ * K * 2 K k £ X ok o c o OI c 0 c X X X XC X X X X Oc O0 X oo 0 X X X OR K KOK K 

REM X Xo oko X oc oc ck £ K Kk KOR c oc X X X X c c X OK kc OK kc ok Xe ox oo xc Xo X X X0 X KOR R Xx X 

PRINT 

REM ******* COMPUTER CONTROLLED MOVEMENT *****x** 

PRINT 

PRINT "THE PRESSURE DATA WILL BE WRITTEN TO FILES ON DRIVE 'A' " 
PRINT 


PRINT "YOU WILL BE ASKED TO INPUT FILE NAMES FOR THESE." 

PRINT 

INPUT. "IS A FORMATTED DISK IN DRIVE 'A'? PRESS ‘ENTER’ TO CONTINUE";D$ 
PRINT 


PRINT 

PRINT 

PRINT " 30k X ko c kc oe ke c c ko kc oko oko c K K KOK KOK kc kc o0 X io 0 ok RK OROR OR R 9 
PRINT " se NOTE ti! KR 
PRINT " ** COMPUTER CONTROLLED MOVEMENT ao 
PRINT " 2. IS PROGRAMMED WITH A ΚΡ 
PRINT " ** DEFAULTED NEGATIVE MOTOR INCREMENT **" 
PRINT " ** (i.e. MOTOR #3 WILL MOVE UPWARD EE 
PRINT " ** BY ENTERING A (+) DISTANCE). ΕΕ 
PRINT " x X K X X X X X X K * kc oko X0 X ok ooo xc xoc kc ook kc c e e e xc ec ee x! 
PRINT 


REM SET INITIAL MOVEMENT DISTANCE AND NUMBER OF DATA POINTS TO ZERO 
HT=0 
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3710 
3720 
3730 
3740 
37509 
3760 
3770 
3780 


WD=0 


DISP 


XPT=0 
YPT=0 
N= 0 

PRINT 
PRINT 
INPUT 


QU GWD.HI 


3790 
3800 
3810 
3820 
3830 
3840 
3850 
3860 
3870 
J880 
3890 
3900 
3910 
3920 
3930 
3940 
3950 
3960 
3970 
39380 
3990 
4000 
4010 
4020 
4030 
4040 
4050 
4060 
4070 
4080 
4090 
4100 
4110 
4120 
4130 
451450 
4150 
4160 
45170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4302 
6304 
4306 
4310 
4320 


PRINT 
INPUT 


YPT=INT(HT /DIST) 


"WHAT IS THE DIMENSION ( X , Y ) (IN INCHES) THAT YOU WANT TO MEASURE 


"WHAT IS THE STEP (IN INCHES) THAT YOU WANT TO MOVE.";DIST 


tal 


XPT=INT(WD /DIST)+ 1 
N-XPT*YPT 


PRINT 
PRINT 


PRINT 


"THERE ARE  ";XPT;" * ";YPT;" - ";N;" POINTS TO BE MEASURED “ 


INPUT "ARE THE NUMBER OF POINTS IS OK.(Y OR N)";C$ 


IF C$="N" 


CLS 

N=XPT 
IF (N 
REM 
B$ 

ES 

X3 

EX$ = 
CLs 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 


ս ի H 


IF LEN(F23) 


CLS 

PRIHT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
IF C$ 
IF C$ 


THEN 3780 


ς 1) OR (N > 99) GOTO 1759 


*** GENERATING STRING STRING SEGMENTS FOR DATA FILE NAMES 

MIDS(STR$(1), 
MIDS(STRS(N), 
"XAXXXX'" 


21 
2): 


REM ** STRING NUMBER "1" 
REM ** ENDING STRING NUMBER "N" 
EDATS 

"DATA FILES WILL BE INCREMENTED FROM:" 


(X$ * B$ * EXSOC "Po ԵԵ Εν LIENS 


"ENTER DATA FILE NAME (6 CHARACTERS MAX -- NO EXTENSION) ";F2$ 


2 6 OR LEN(F2$) < 1 GOTO 4030 


N; "DATA FILES WILL BE GENERATED AND INCREMENTED‘ AS FOLLOWS?” 


(Բ25 + BSet EXS ) 5 (To i (B2t + 1 Խա Է) 


"ARE THE NUMBER OF POINTS AND FILE NAMES OK.(Y OR N)"; C$ 
= "N" GOTO 3780 | 
= "y" GOTO 4190 


GOTO 4140 


CLS 
PRINT 
PRINT 


REM SET INITIAL 


POSITION DATA 


X(1)=-DIST 
Y(1)=-DIST 
FOR IX=2 TO XPT*1 


τν. 


0 


NEXT IX 
FOR —2 15 ΡΠ 


Y(JY)s 


0 


NEXT JY 


FOR I= 


τι 
12 O 
1 -ն 


FOR J= 
MOTOR CP-CONTROLLED MOTOR MOVEMENT 


REM 


1 TO XPT 


ο τει 
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4330 11=0 

5340 I2-0 

4350 13-0 

4360 REM EACH POINT TAKE 10 TIMES READINGS 

&370 X(I+1)=X(I)+DIST 

4380 XPT(J)=X(I+1) 

4390 Y(J+1)=Y(J)+DIST 

4400 YPT(J)=Y(J+1) E 
4405 INPUT " ADJUST THE WHEEL TO MAKE THE P2 =P3,INPUT THE YAW ANGLE"; YAW(J) 
4408 PRINT 

4410 INPUT " PRESS 'ENTER' TO START THE MEASUREMENT" ;MOVE$ 
4420 REM 

4430 REM READ FIVE CHANNELS AND DISPLAY THE DATA 

4440 REM 

4450 STEPPER=4 

4460 SWITCH = 3 

4470 HOMER=8 
4480 DELAY1 = .1 

4490 DELAY2 = 1 

4500 REM SET THE S.V PORT TO էգ 

4510 FOR IL=1 TO 3 

4520 THYME = TIMER 

4530 CALL OUTPUT(RELAY.ACT.O1,STEPPER) 

4540 CHKTIME = TIMER 

4550 IF CHKTIME < (THYME + DELAY1) GOTO 4540 

4560 CALL OPEN. CHANNEL(RELAY.ACT.01,SWITCH) 

4570 CLS 

4580 NEXT IL 

4590 PRINT 

4600 PRINT " NOW IS POINT ";J 

4610 REM START MEASURE FROM PORT & TO PORT 8 

4620 FOR JJ=1 TO 5 

4630 CALL OUTPUT(RELAY.ACT.01,STEPPER) 

4640 CHKTIME = TIMER ` 
5650 IF CHKTIME « (THYME * DELAY2) GOTO 4640 

4660 REM EACH PORT SAMPLE 10 TIMES 

4670 FOR II=1 TO 10 . 

4680 ROUT=1 

4690 CALL OUTPUT(RELAY.MUX.01,ROUT) 

4700 CALL MEASURE(DMM.01,VOLTS) յ 
710 PA(II,JJ)=VOLTS 

4720 NEXT II | 
4730 CALL OPEN.CHANNEL(RELAY.ACT.01,SWITCH) 

4740 IF JJ=5 THEN 4760 

4750 NEXT JJ 

4760 REM HOME THE S.V PORT TO #48 

4770 CALL OUTPUT(RELAY.ACT.01,HOMER) 

4780 CALL OPEN. CHANNEL(RELAY.ACT.01,HOMER ) 

4790 REM 

4800 REM DISPLAY THE SAMPLE DATA 

4810 REM 

4820 PRINT HEAD1$ 

4830 FOR IS= 1 TO 10 

4840 PRINT USING FORMAT$;1S,XPT(J),YPT(J),PA(IS,1),PA(IS,2),PA(IS,3),PA(IS,4&),PA 
(IS,5),YAW(J) 

4850 NEXT IS 

4860 REM 

4870 REM AVERAGE THE DATA 

4880 REM 

4890 FOR JA = 1 TO 5 

4900 TOTAL = 0 

4910 FOR IA = 1 TO 10 

GO20, TOTAL = TOTAL t PACIA,JA) 

4930 NEXT IA : 
4940 AVERAGE = TOTAL 710 

4950 P(JA)=AVERAGE 

4960 NEXT JA 


Սն 


4970 PRINT 

4980 PRINT “THE AVERAGES ARE: " 

5000 PRINT HEAD1$ 

5010 FOR JD=1 TO 5 

5020 PP(J,JD)=P(JD) 

5030 NEXT JD 

S040 PRINT USING FORMATS$;J,XPT(J), YPT(J) , PP(J, 3), PP(J, 2), PP(J, 3), PP(J, 4), PP(J, 5) 
ΥΛΗ (0) 

5045 PRINT 

5048 PRINT USING "THE NULLING ERROR IS +#.88#48";PP(J,3)-PP(J,2) 

5049 PRINT 

5050 PRINT "DO YOU WANT RE-MEASURE AGAIN (Y / N)" 

5060 PRINT 

5062 PRINT "IF 'Y' WILL RE-SAMPLE AGAIN." 

5064 PRINT 

5070 INPUT "IF 'N' WILL MOVE THE TRAVERSER STEP UPWARD (WAIT 7 SEC )";C3 
5075 PRINT 

5080 IF C$-"Y" THEN 4405 

5082 IF C$="N" THEN 5090 

508% GO TO 5070 

5090 IF J-YPT THEN 5160 

5100 REM 

5110 REM MOVE THE TRAVERSER STEP UPWARD. 

5120 REM 

5130 13=-DIST 

5140 GOSUB 3410 

5150 NEXT J 

5160 REM*** STORE DATA BEFORE NEXT SAMPLE*** 

5170 OPEN "A:NRAW.DAT" FOR OUTPUT AS #2 

5180 PRINT #2 ,HEADI$ 1 
5190 FOR ID-1 TO YPT 

5200 PRINT 42 ,USING FORMAT$;ID,XPT(ID),YPT( ID), PP(ID, 1), PPCID,2) , PPC ID,3) , PP(ID 
4), PPCID, 5), YAWCID) 

5210 NEXT ID 

5220 CLOSE #2 
5230 REM *** GENERATING INCREMENTED DATA FILE NAME 

5240 IF (I > 10) OR (I = 10) THEN I$ = MID$(STRS(I), 2) 

5250 IF (I « 10) THEN I$ = (MIDS(STR$(0), 2) + MIDS(STRS(I), 2)) | 
5260 FI2$ - (F2$ + I$ + EX$) 

5270 PRINT | 

5280 PRINT " WRITING DATA FILE "; FI2$ 

5290 DF2$-RN2$4FI2$ 

5300 REM ** RENAME DATA FILE 

5310 SHELL DF2$ 

5320 REM | 

5330 REM MOVE THE TRAVERSER TO THE NEXT SAMPLE POSITION 

5340 REM 

5350 PRINT 

5360 IF I-XPT THEN 5430 

5370 INPUT "THEN PRESS 'ENTER' FOR NEXT COLUMN SAMPLE( 90 SEC) ";MOVE$ 
5390 I2--DIST 

5400 I13=HT 

5410 GOSUB 3410 

5420 NEXT I 

5430 CLS 

5440 PRINT "ALL MOVEMENTS COMPLETE" 

5450 PRINT 

5460 PRINT 

5470 PRINT "YOU WANT TO REPOSITION TRAVERSER FOR ANOTHER MOVEMENT (Y OR N)?" 
5480 PRINT 

5490 PRINT "IF 'Y', THE PROGRAM WILL TAKE YOU TO MANUAL CONTROL SUBROUTINE." 
5500 PRINT "IF 'N', THE PROGRAM WILL END," 

5510 PRINT 

5520 INPUT "ANOTHER MOVEMENT" ;R$ 

5530 IF R$ = "Y" THEN 1370 

5540 IF R$ = "N" THEN 3370 
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APPENDIX B. CALP PROGRAM 


DEF SEG:CLEAR ,&HFEOO:GOTO 4 'Begin PCIB Program Shell 
GOTO 1000 ' User program 

GOTO 900 ' Error handling 

I-&HFEOO ' Copyright Hewlett-Packard 19845,1985 
PCIB.DIR$-ENVIRONS$("PCIB") 

IP=PCIB.DIRE+"\PCIBILC. BLD” 

BLOAD 18,1 

CALL I(PCIB.DIR# ,1%,J3J%):PCIB.SEG=1% 

IF J%=0 THEN GOTO 13 


PRINT "Unable to load."; 
PRINT " (Error 94 da) 


END 
t 


DEF SEG=PCIB.SEG:O.S=5:C.S=10:I,V=15 
MG z20: Gi Ps25 DETE 0 


GET .MEM=35:L.5=40: PANELS=45: DEF .ERR=50 

PCIB. ERR#=STRINGS (64,32) : PCIB.NAMES=STRINGS$(16, 32) 

CALL DEF. ERR(PCIB. ERR, PCIB. ERRS, PCIB. NAMES. PCIU.GLBERR) : PCIR.BASERR=255 
ON ERROR GOTO 3 

J=- 


I$sePCIBDIP$SS"NPCIB.SYN" 

CAELO SOLES) 

IF PCIB.ERR<>20 THEN ERRORWPGTB.BASERR 

1Ξ0 

CHDESPSVUOISREADZREGISTEBRMBEAD.SEDEIDS.DEEINE,TNITIALIZE.SYSTEM) 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CAEL I1.V( (I BNABLE. SYSTEM (DISCKHGE.SsYSTEM, INITIALIZE, POWER.ON) 

IF PCIBL ERR<>0 THEN ERROR PCIBBASERR 

CALL IS V(I MEASURE. OUTPUT STARI  NALT) 

IF PCIÐB.ERR<>0 THEN ERROR PCIB I PASEPR ! 
CALL I.V(I,ENABLE.INT.TRIGGER,DISABLE.INT.TRIGGER, ENABLB. OUTPUT , DISABLE.OUTPU 


LECPCLIBOERRCGODSPHENCERBRORSPGEB-HASERR 

CALE Sra CHECK. DONE (CET. StAtUs, SET “FUNCTION, SET .RANGE) 
IF PCIBSERR<>O0 THEN ERROR PCIB BASERR 

CALL ΕΕ ο SEP -MODE ΠΕΙΤΕ GAL READ. CAL, STORE.CAL) 
IF PCIÐ.ERR<>0 THEN ERRORTPCIB.BASERR 

CALL I. VCI DELAY; SAVE: SISTEM JJ? 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

TEJ 

CARLE 1I: VCI, SET.:GNTETIME SET. SpHPLES SET.SLOPF , SET.SOURCE) 
IP POTI ΕΕ SOT THEN ERROR ἘΕΙΡΕΠΞΕΓΕ 

ΠΤΙ T EFREOUENCY AUTO. FREQ, PERIOD, AUTO. PER) 
IF PCIB.ERR<20 THEN ERROR PCIBÐB.BASERR 

GALL T.C(Ul, INTERVAL, RATIO, FOTALIZE,RIOOMILLI ) 

IF PCIB.ERR<>0O THEN ERROR, PCIB.HASERR 
ΠΙΟ... Ισ RIKILO) 

DF ΡΟΙΕ ΕΡΕ ο ΤΕΝ ERROR PCIB BASERR 

CALL I.C(1I,RIOMEGA,R1OOMEGA, CHAN.A,CHAN.B) 

IF PCIB.ERR<>0 THEN ERROR PCIB. BASERR 

CALL 17CU1l,POSITIVFE..NEGATIVE, COMN, SEPARATE) 

DE Pel BR OEBRRCOOSTHEHSERPSROPGTBOBASERE 

152 

153 
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106 
107 
108 
109 


CALL I- V(I, ZERO OMS SELT: S ERD J «I 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 

CALL I.C(I,DCVOLBIS.ACVOLTS,OHMS,R200HMDISELR) 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL I.C(I,R2,R20,R200,R2KILO) 

IF PCIB.ERRC^D THEN ERROR PCIE. BASERR 

CALL I.C(I,R20KILO,R200KILO,R2MEGA,R20MEGA ) 

IF PCIB.ERR«X»0 THEN ERROR PCIB.BASERR 

CALL I.C(I,AUTOAW R2 S PB12. 5 p) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 

1-4 

CALL I.V(1,SET.COMPLEMENT ,SET.DRIVER,OUTPUT.NO. WAIT, ENABLE SHARDS el 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL I.V(1,DISABLE.HANDSHAKE, SFT.THRESHOLD, SET .START. BIT, SE] Πα... 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL I. V(I, SET: LOGIC., SENSE; J JA 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL I.C(J,POSITIVE,NEGATIVE,TWOS,UNSIGNED) 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL σος το πο πι 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 

CALE I.CULVRZSBE3S IO RO 

IF PCIB.ERR<>0 THEN ERROR PC] B{BASERR 

Kn Kee 

IF PCJ] B.ERR<>0 THEN ERROR PCIB. BASERR 

ΕΕΒΕ Ίο ΕΙ. ΕΙ οι P 13) 

IF PCIB.ERR<>0 THEN ERPOR PCIB.BASERR 

CALL. EES EH 

IF PCIB.ERR<>O THEN ERROR PCLBOBASERR 

156 

CALL 1.V(.1,SET.FREQUENCY SEP AMPLITUDE SET OF FSER SET oS 7 (ite 

JF PCIB.ERPR<50 THEN ERRORSPCIE, BNOeRR 

CALL I.:V(1,SET. BURST {Coun Jn) D) 

IF PCIB.ERR«»0 THEN ERROR FCIB.BASERR 

CALL 1.C(J],SINE,SQUARE,. TRIANGLE , CONTINUOUS ) ' 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 

CALL -L.C GATED, BURST AJ 

IF PCIB.ERR<>0 THEN ERROR PCIB T BASERR 

157 

CALL I.V(1,AUTOSCALE, CALIBRATE .SET SEMSITIVITY SET ΕΕ OPP Ps?) 

IF PCIB.ERR<>0 THEN ERROBRSEPCIBOBASERI 

CALL I.V(I,SET.COUPLING,SET.POCARTIY SET SWEEPS Ρε τι 

IF PCIB.ERR“<>O THEN ERROR PCIB.BASERR 

CALL 1.V(1,SET.TRIG. SOURCE,SET.TRIG.SLOPE,SET.TRIGJLEVEGUSGET AJ Ra i Մ 
IF PCIB.-ERR<>0 THEN ERROR PCIÐB BASERR 

CALL 1.V(1,GET.SINGLE.UF,GET.TUO.VE,GET.VERT INPO, GET TIME RAShHI OJ 
LE PCIB.ERR<>0 THEN ERROR PCIE BASERR 

CALL I.V(I,GET.TRJG.INFO,CALC.WFVOLT,CALC.WFTIME,CALC.WF.STATS) 
IF PCIB.ERR<>0 THEN ERROR PC IB OPP ASE 

CALL I.V(I,CALC.RISETIMFE,CALC.FALLTIMI,CALC.PERIOD ,CALC. FREQUENCY) 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL I.V(I,CALC.PLUSWIDTH, CALC. MINUSWIDTH, CALC. OVERSHOOT , CALC . PRESHOOT) 
TP PCIB.ERR<>0 THEN ERROR PCIBTBASERRK 

CALL I.V(I,CALC.PK.TO.PK,SET.TIMEOUT,SCOPE.START,MEASURE.SINGLE.WF ) 
IF PCIB.ERR<>0 THEN ERRORCPCLIBEBASERR 

CALL I.V(I,MEASURE.TWO.WF,J,J,J) 
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110 
111 
112 
JJ) 3 
114% 
115 
LEG 
1117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
ΤΠ 
134 
1535 
136 
1317 
138 
139 
140 
141 
142 
ΠΗ 5 
144% 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
219 
156 
157 
158 
159 
160 
161 
162 
163 
165 
800 
801 
899 


IF PCIB.ERR«^O THEN ERROR PCIB.BASERR 

CALL I.C(I,R310HANO,R100NANO,RI1MICRO,R10MI CRO) 
IF PCIB.ERF«^O0 THEN ERROR PCIB.BASERR 

CAEDE I.CUITRIOOMICRORIMIEBEI.RIOMIELI,RIOOMILLEI ) 
IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 

CALL I.C(I,R1,R10,R20NHANO,R200NANO) 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

CALL I.C(I,R2MICRO,BR20MICRO,R200MICRO, R2MILLI ) 
IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 

CALL I.CCI.R20MPEEI.R200MIBEI,.R2;R20) 

IF PCIB.ERRKSO THEN ERROR PCIB.BASERR 

CALL I.C(1,RSONANO,RSOONANO,RSMICRO,RSOMICRO) 
IF PCIB.ERR«»0 THEN ERROR PCIB.PASERR 

CALL I.C(ISR500MICBRBOSRS5MIBET,RSOMILET.RSOUMILLI) 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 

Շու Ἱ οι R5,RBR50.CIIRH IA CHAN. B) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 

CALL I.C(I, EXTERHAL,POSITIVE,NEGATIVE, RC) 

IF PCIB.ERR<>O THEN ERROR PC1B.BASERR 

CALL Ἱ ο EEH 
IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 

CALL I.C(I,X1,X10,STAHDARD,AVEFAGE ) 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

I-8 

CALL 1:V (1. OPEN CHANNEL. CLOSE. CIINMMEL, J..)) 

IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 

CALE C.S 

IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
IT-PCIB.DIRT*" NPCIB.PLD" 

CALL L.P (IS) 

IF PCIB.ERR«^O THEN ERPOR PCIB.BASERR 
l¢="DMM.01":I=3:J=0:K=0:L=1 

CALL DEFINE(DM1.01,I$,I,J,K,L) 

IF PCIB.ERP<>0 THEN ERROR PCIB.BASERR 
I$="Func.Gen.01":I=6:J=0:K=1:L=1 

CALL DEFINECFUIIC.GEN.0O1, Me, 1.J,K,L) 

IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
I$-"Scope.01":1-27:J-7-0:K-2:L71 

CALCE DEFINE (SCOPE Or TIS T J K EC) 

IF PCIB.ERR<>0 THEN ERROR PCIR.BASERR 
DLt-"Counter.01":7-1:J-0:K*53:L*1 

ՇՈՒՆ DEFIHE(COUNTER.01,I$,I,J,K,L) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 

Ἱ DIPINJOI JEN IO Ἱ-] 

CREL ΡΕ ος ο. οἳ IS ο ο D) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 
I$="Dig.Out.01":I=4:J=1:K=4:L=1 

CREL DEFINE(DIC OUT Or DET JALKE) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 
I$="Relay.Act.01":I=8:J=0:K=5:L=1 

CALL DEFINE(RELAY.ACT:01,r$,I.J,K,L) 

IF PCIB.ERR«»0 THEN ERROR PCIB.BASERR 
I$f-"Relay.Mux.01":1-7-2:J720:K726:L71 

CALL DEEIBHMECBELAY O MUK.O1 1 L.J K,;Db) 

IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 
IS=ENVIRONE ( "PANELS" )+"\PANELS. EXE" 

CIS S IS) 

GOTO 2 ç 
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900 IF ERR=PCIB.BASERR THEN GOTO 903 

901 PRINT "BASIC error $";ERR;" occurred in line ";ERL 

902 STOP 

903 TMPFERR=PCIB.ERR: IF TMPERR=0 THEN TMPERR-PCIB.GLBERR 

904 PRINT "PC Instrument error *";TMPERR;" detected at line ";ERL 

905 PRINT “Error: ";PCIB.EHRS 

906 IF LEFT$(PCIB.NAMES$,1)«»CHR$(32) THEH PRINT "Instrument: ";PCIB.NAMES$ 

907 STOP 

908 COMMON PCIB.DIR$,PCIB.SEG 

909 COMMON LD.FILE,GET.MEM, PANELS, DEF.ERR 

910 COMMON PCIB.BASERER,PCIB.ERR,PCIB.ERR$,PCIB.NAME£$,PCIB.GLBERR 

911 COMMON READ.REGISTER,READ.SELFID,DEFINE.INITIALIZE.SYSTEM, ENABLE.SYSTEM,DISA 
BLE.SYSTEM,IHITIALIZE,POWER.OH,MEASUPRE,OUTPUT,START,HALT,ENABLE.INT.TRIGGER,DISA 
BLE.INT.TRIGGER,ENABLE.OUTPUT,DISABLE.OUTPUT,CHECK.DONE,GET.STATUS 

912 COMMON SET.FUNCTION,SET.RANGE,SET.MODE,WRITE.CAI.,READ.CAL,STORE.CAL,DELAY,SA 
VE.SYSTEM,SET.GATETIME,SET. SAMPLES, SET. SLOPE,SET. SOURCE, ZERO. OHMS SET SEE r 
COMPLEMENT,SET.DRIVER,OUTPUT.HO.WAIT,EHABLE.HANDSHAKE,DISABLE.HANDSHAKE 

913 COMMON SET.THRESHOLD,SET.START.BIT,SET.NUM.BITS,SET.LOGIC.SENSE,.,SETOBREGDUBIE 
Y,SET.AMPLITUDE.SET.OFFSET,SE!I.SYMMETRY,SET.BURST.COUNT , AUTOSCALE, CALIBRATE a 
SENSITIVITY,SET.VERT.OFFSET, SET.COUPLING, SET. POLARITY ,SET.SWEEPSPEED 

914 COMMON SET.DELAY,SET.TRIG.SOURCE,SET.TRIG.SLOPE,SET.TRIG.LEVEL,SET.TRIG.MODE 
,GET.SINGLE.WF,GET.TWO.WF,GET.VERT.INFO, GET. TIMEBASE. INFO, GET. TRIG. IHFOS; CAPERE 
OLT,CALC.WFTIME,CALC.WE-STAIS,CALC.RISETIME,CADCOEALLAMME CALC TERI OPE 

915 COMMON CALC.FREQUENCY,CALC.PLUSWIDTH,CALC.MINUSWIDTH,CALC.OVERSHOOT, CALC.PRE 
SHOOT,CALC.PK.TO.PK,SET.TIMEOUT,SCOPE.START.MEASURE.SINGLE.WF,MEASURE-TWODWENMDEE 
H.CHAMMEL,CLOSE.CHBANMEI 

916 COMMON FREQUENCY,AUTO.FREO,PERIOD.,AUTO.PER,INTERVAL,RATIO, TOTALIZE,R100MILLI 
,R1.R10,R100,R1FILO,R10MEGA,R100MEGA,CHAN.A,CHAH.B, POSITIVE, NEGATIVE, COMN, SEPARA 
TE,DCVOLTS,ACVOLTS,OHMS.R200MIDEI.,R2.R20.R200,R2KTIDELOSR2OPRUEO.R200PIEU 

917 COMMON R2MEGA,R20MEGA,AUTOM,R2.5,R12.5,POSITIVE, NEGATIVE, TWOS,UNSIGÜHED,OC, TT 
L,RO,R1,.R2,R3,Ré&,RS,R6,R7,RB,RI,R1O,R11,R1i2,R13,PR14,RiS5,R1i6,SINE, SQUARE, TRAV eee 
»CONTINUOUS, GATED, BURST, R1IONANO, RIQOONANO, RIMICRO,PiOMICRO,RiIOQOMICRO 

918 COMMON RIMILLI,RiIOMILLI,RIOOMILLI,R1,R10,R20NANO, R2ZQ00NANO, R2MICRO,R2Z0MICRO,R 
200MICRO,R2MILLI,R2Z0MILLI,R2Z00MILLI,R2,R20, RSONANO, RSOONANO, RSMICRO, RSOMICRO,R50O 
OMICRO,R5MILLI,RSOMILLI,RSOOMILLI,R5,R50, CHAN.A, CHAN.B, EXTERNAL, POSITIVE 

913 COMMON NEGATIVE,AC.DC,TRIGGERED,AUIO.TPIG,AUTO. LEVEL, X1,%10, STANDARD, AVERAGE 
920 COMMON DMM.01,FUNC.GEN.01,SCOPE.01,COUNTER.0O1,DIG.IN.01,DIG.OUT.01,REDAY AG. 
- 01, BELAY 7HMUx 701 

999 'End PCIB Program Shell 

1000 REM 

1010 REM This step initialzes the HP system 

10οέο τς 

1030 OPTION BASE 1 

1040 DIM P(10),PA(50,6) ,PP(50.6), 114401 Նոն(ն0) 

1050 CALL IHITIALIZE.SYSTEM(PGMSHEL.HPC) 

1060 REM 

1070 REM All PC devices now have an initial state 

1080 REM Set function on the DMM and Relay MUX 

1090 REM 

1100 CALL SET.FUNCTION(DMM.01,DCVOLTS) 

1110 CALL SET. RANGE(DMM.01,AUTOM) 

1120 CALL DISABLE. PNT. TRIGGER(DRM: 01) 

1130 CALL ENABLE.OUTPUT(RELAY.MUX.01) 

1140 FORMATS="#8 {1.185} FARRAR 15.111} OA RH OR OR 

1200 FOR I*1 TO 10 

1210 CAL(1I)=0.0 
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1220 NEXT I 

1510 REM 

1520 REM READ THE VOLTAGE OF 48TH CHANNEL AND DISPLAY THE DATA 
1530 REM 

1540 PRINT " CHOOSE 6 POINTS" 

1550 PRINT 

1550 PRINT "THE CALIBRATION WILL BE STORES IN 'CAL.DAT'" 
1560 REM 

1570 REM Begin sampling loop 

1580 REM 

1600 FOR J=1 TO 1 

1610 PRINT 

1630 FOR JJ=1 TO 6 

1631 INPUT "INPUT THE CALIBRATION PRESSURE";CAL(JJ ) 

1632 INPUT "PRESS 'ENTER' TO START MEASUREMENT";MOVE$ 
1640 FOR JI=1 TO 10 

1650 ROUT=1 

1660 CALL OUTPUT(RELAY.MUX.01,ROUT) 

1670 CALL MEASURE (DMM.01,VOLTS) 

1680 PA(II,JJ)-VOLTS 

1690 NEXT II 

1700 IF JJ=6 THEN 1740 

1730 NEXT JJ 

1740 REM 

1750 REM DISPLAY THE SAMPLE DATA 

1760 REM 

1780 FOR IS* 1 TO 10 

1790 PRINT USING FORMAT$;IS,PA(1S,1)2,PA(IS,2),PA(IS,3),PA( IS, 65), PAC(IS, 5), PAC(IS,6 


16800 NEXT IS 

1910 REM 

1820 REM AVERAGE THE DATA | 
1630 REM 

1840 FOR JA = 1 TO 6 

r850 TOTAL = 0 ' 
1860 FOR IA = 1 TO 10 

1870 TOTAL - TOTAL * PA(IA,JA) 

1880 NEXT IA 

1890 AVERAGE - TOTAL /10 

1900 P(JA)-AVERAGE 

1920 NEXT JA 

1930 PRINT 

1960 PRINT "THE AVERAGE ARE: " 

2000 FOR JD=1 TO 6 

2010 PE(J,JD)=P(JD) 

2020 NEXT JD 

55 PRINT USING FORMATS :J3,PPCJ,1),FP(J,2),PP(J.3),PP(0,4),PP(J,5),PP(0,6) 
2070 PRINT 

2080 INPUT "DO YOU WANT RE-MEASURE AGAIN ? (Y / N)";C$ 

2090 IF. C$-"Y" THEN. 1580 

2101 REM*** STORE DATA BEFORE NEXT SAMPLE*** 

2102 OPEN "A:\CAL.DAT" FOR OUTPUT AS #2 

2106 FOR ID=1 TO 6 

2107 PRINT 82,USING FORMAT$;ID,PP(J,ID),CAL(ID) 

2108 NEXT ID 

2709 CLOSE 12 

2210 NEXT J 
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APPENDIX C. CONVERT PROGRAM 


X X k 3k k K OR 0k oc kc X c kc X0 kk OX ok kc kc kc kc ok X X kc oc kc Xe X kc OOOO oto 5R OR OR OR OR ËR OR OR OR OR OR R R R k A 9 8 9 tote st 
* THIS PROGRAM CONVERTS THE VOLTAGE OF TRANSDUCER INTO PHYSICAL * 


* PRESSURE,VELOCITY,YAW ANGLE AND PITCH ANGLE. THOSE DATA ARE 


* USED FOR PLOT PROGRAM LATER. 
* X X Xk X XK kk k 3k k 3k OR Y OR OR OR OR R Y OR OR R OROR OR OR R R ORO R OR R KOR OR R OR ORO OR R R R R R ROF k R R R R R R KOR R R R R k k 


* 


* X *%* % H 


CHARACTER*12 FNAME 
CHARACTER*12 NAME 
CHARACTER*2 A(50) 
CHARACTER*80 ST 

REAL K,INTR 

INTEGER COLS,RWS,DTPTS 


DATA A/'01','02','03' , "On MINUS ան Սա աա ի 
10, 11,127. 13221015. Ան ի 
19; 20’, 121,22!» ΤΟΝ ο οσο. 
"28°" 591, 5 Ὁ. ՍԱ ն 
3738 "3097 ΠΠ’, 4] Ὁ ο ου ն 
երո Կ կն նն 2 րԱ 

WRITE (*,'(A\)') ' # OF COLS (AWAY FROM MSL) = ' 

READ <(*).' (15) COLS 

WRITE (*,'(A\)') ᾿ # OF DATA PTS IN A COL (UP/DOWN) 

READ (*,'(15)5) RWS 

WRITE (*,'(AX)') ' DATA FILE NAME? (IE ROOJAZXX. DAT) 

READ): ԱՆ 

WRITE 417. 1 Ա τι εὖ ο) i 

READ (Էն 9-1 

WRITE 17, CANON Բր. πα 2] 

READ ει οἱ PE 

WRITE ("6.27 Խան Ք Ս | 

READ (ET 

ՍՔ 11, (2. (Ս ΤΕΕ 

րոն (Es. ՆԲ 

ΜΕΗΙΙΕ ΙΤ ΠΤ, ἨΜΙΡΕ νυν, 

REND lee 

WRITE (*,' (AN)') © SLOPE FOR DELTAP S(T? 51 C E 

READ ΙΙ; (REGJ) SONE 

WRITE (*,'(AX)') ' INTERCEPT FOR DELTAP (F9.6) s ' 

READ (32, '(F9.6)') INIR 

WRITE (*,'(NN)') ' QM1 FACTOR (F4.2) = ' 

READ (* '(Fh&,2)') QMHIEñnG 

WRITE ere 

READ (*,'(F5.,2)') XOFE 

WRITE (*,'(AN) y) 59.1. 

ΒΕΑΛΡ ΟΕ (CES. 2) 1] vOEE 


PATM-(PI*PF)*35.363] 
R=1716.5 

E=0.0123 
T=(TI+TF)/2. +460 
RO-PATM/(R*T) 
DTPTS-RWS*COLS 
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CONVERT THE PRESSURE UNIT FROM inHg TO psf 


' 


* 
* 


* OPEN A NEW FILE TO STORE THE REDUCED DATA 
OPEN(2.FTEE-'"RESUDT.DAT',STATUS-'OLD') 
WRIIE(2,222) DTPTS 

2 2 2 FORMAT (15) 

* OPEN A SEQUENTIAL OF DATA FILE 

DO 2071 1 CODS 
NAME(7:8)=A(T) 
FNAME=NAME 
OPEN( 1, FILE-FNAME) 

READ(1,100,END=20)ST 


100 FORMAT(A65 ) 
15 READ(1,1000, END 30)NO, X, Y, V1,V2,V3,V4, V5, BETA 
1000 FORMAT(I27FEF7 27;F6%°2 7 5EFE9 3 7 FEF8.2) 


* CONVERT THE VOLTAGE TO PRESSURE IN LBF/FT**2 
P1-DELTAP(V1,SLOPE,INTR)*2.0475*PATM 
P2-DELTAP(V2,SLOPE, INTR)*2.0475*PATM 
P3-DELTAP(V3,SLOPE,INTR)*2.04575*PATM 
P4=DELTAP(V4,SLOPE, INTR)*2.0475+PATM 
P5-DELTAP(V5,SLOPE, INTR)*2.0475*PATM 

* CALCULATE THE PITCH ANGLE IN DEGREES 
P=(P4-P5)/(P1-P2) 

ALPHA=FPITCH(P) 

* CALCULATE THE VELOCITY IN FT/SEC 
YSLOP=FYSLOP(ALPHA ) 
VELM-SQRT((2*YSLOP*(P1-P2))/(RO*K)) 
VEL-VELM*(1*E) 

* CALCULATE THE LOCAL DYNAMIC PRESSURE 
OM1=OM1FAC*2.0475/K 
QM-RO*VEL**2/2. 

Q1-QM1*(1*2*E) 
Q-QM*(1*42*E) 

* CALCULATE THE YAW ANGLE IN DEGREES 
YAW-FYAW(BETA*5.0) 

* CALCULATE THE VELOCITY COMPONENTS 
BETAR-YAW*.017453 
ALPHAR-(ALPHA-17.942)*.017453 
VELY=VEL*SIN(ALPHAR) 
VELX=VEL*COS(ALPHAR)*SIN(BETAR) 

* CALCULATE THE TOTAL PRESSURE IN LBF/IN**2 
PTC=FPT (ALPHA) 

PT1=P1-0'%PTC 
ENEE u 
CPT-(PT1-PATM-Q1)/01 

* CALCULATE THE STATIC PRESSURE IN LBF/IN**2 
բտ1-Ք11-0 
PS=PS1/16% 

CPS-(PS1-PATM)/Q13 

* WRITE ZERO VALUES IF VELOCITY IS TOO HIGH (BAD) 

| IF(VEL.GT.200.0) THEN 


670.0 

ΤΕ O000)=- X πο ԵՒՆ ο ο ὁ ασ .,ᾱ 

ELSE 

WRITECZ p20U00 ) -X*XOPRF .Y*YORF (VEL. VELX, VELY, YAW, 
C ALPHA-17.942,PT,CPT, PS, CPS 
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2000 FORMAT(11F10.3) 
ENDIF 
GO TO 15 
30 CLOSE(1) 
20 CONTINUE 
CLOSE(2) 
STOP 
END 
KKK KK KKK KKK TKK RK KKK RRR RK KKK KEK RK KKK KKK RK KERR RK KKK kt KKK KKK KK AK 
* THIS FUNCTION CONVERTS THE VOLTAGE TO PHYSICAL PRESSURE 
FUNCTION DELTAP(X,SLOPE, INTR) 
REAL INTR 
DELTAP=X*SLOPE+INTR 
END 
KKK KTR KK KR RK RK KKK KKK ERK KE KKK KEE RR KKK KT KT KKK KKK KKK KKK KKK KY F 
* THIS FUNCTION CALCULATES THE PITCH ANGLE 
FUNCTION FPITCH(X) 
FPITCH=3.759+53.7568*X-1.3085*X**2-1.6583*X**3 
* -0.8061*X**4+16.5115*X**5 
END 
2szs2z29229o922222222229222292222222222222922229 22222292 2222222225 22:222252222 
* THIS FUNCTION CALCULATES THE VELOCITY PRESSURE COEFFICENT 
FUNCTION FYSLOP(X) 
IF(X.LT.-10)THEN 
FYSLOP-0.981-0.0102*X-3.000E-4*X* *2-2.500E- 6*X* *3 
ELSE IF((X.GE. -10) . AND. (X. LE. 10) ) THEN 
FYSLOP-0.98-0.006*X*2.000E-4*X**2 
ELSE 
FYSLOP-0.9801-0.0035*X-1.1503E-4*X**245.833E- 6*X**3 
END IF 
END 
* oo ox o x o xx x KEKE KKK EE KKK KEKE ER EK KKK EEK EKA KKK KKK KS 
* THIS FUNCTION CALCULATES THE YAW ANGLE 
FUNCTION FYAW(X) i 
IF((X.GE.0).AND.(X.LE.180)) THEN 


FYAW=-X 
ELSE 
FYAW=360-X 
END IF 
END 


ΚΚ ΧΧΧ ΧΑ κα ἐ ΧΧΧ ΧΧΧ ΧΧΧ 1 ERR KEKE KKK Ke 


* THIS FUNCTION CALCULATES THE TOTAL PRESSURE SGerrrICrenNt 
FUNCTION FPT(X) 
PP IGE sd © Peery 
FPT=-0.01 
ELSE IF((X%.GI.-30), AND 2 τε ο THEN 
i FPT-0.02+1,00R +372) 
ELSE IF((X,GE. 520) ANDACX DE JONET JEN 
EDI 
ELSE 
FPTs0.03-1.00E-3*X 
END IF 
END 
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APPENDIX D. RESULT 00.DAT 


84 
1 .000 000 122209 158 18.214 155002 2017 
.875 = 060 
2 .000 900 1220155 * 72:500 18.388 2590072 5 8/2 
2875 ΟΕ 
3 000 1-060 13 20207 = 000 18.254 159002 2017 
2075 soe 
4 .000 1.500 122-99 =7 2000 15 297 15.002 .016 
2875 «061 
S 2000 2.000 1227310 =7 000 19.433 15002 TOTE 
875 ΠΕΙ 
6 2000 2250 0 1220155 τυ 18.388 T5 002 .016 
875 ΠΕΙ 
7 «000 32.000 1225293 “15000 18:302 1°. 092 018 
875 DS 
8 000 3,300 1220359 =7 2.000 181.287 155002 “018 
875 2.0672 
3 .000 4.000 122.287 «ԹԱՍ 18:357 1 ου 017 
875 “2062 
TO .000 &.500 129595871 =5 r000 1755651 15:002 . 017 
868 ΞΘ 
11 -000 57 000 12599793 -5:000 19457373 15 002 «03:8 
868 EU 
12 ους 92500 12598987 τυ 17.614 15.002 . 019 
867 x» 
1 2500 .000 12505938 5952000 17.784 JS. 002 017 
867 e 
2 «ՏԱՍ z300 12557323 “ 5000 1285723 15:002 018 
868 "0122 i 
3 SHU 1.000 126131] =>. 900 17.814 125.002 -015 
866 πο 
& ΤΟ 1.500 12511113 =5 000 17.714 15002 “UJE 
867 El 
5 -500 2000 1262015 25:901 17.764 19.002 015 
867 SIE 
6 . 500 2500 1257812 25000 17.655 15-002 2013 
867 - 1285 
7 «ՅԱՍ οσο 125.624 ms 98 17.669 15.001 2012 
867 [2426 
8 200 359500 IZ m1 5.000 125570 15.002 .015 
867 = 126 
9 -500 4.000 1259777 "ο 68 17.699 15.002 -014 
867 ο αὶ 
10 «5900 &.500 125.944 -5.000 17.601 τς 002 ΠΠ 
867 Ig 
11 2300 5:000 126.524 -5.000 17.454 15.002 015 
2965 -.140 
12 2500 5.500 1259918 «5:900 17.784 15.002 οι” 
«86 τα. 
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«866 


867 


.866 


065 


.866 


.866 


| ο ο ο ՐԾ 1 00 4 P9 JI | I FE O04 EB ԵԹ ԻԹ ԻԹ ԻԹ ԵԹԹ Թ Ո Թ Ե Թ Ա Թ յ թ յ Թ IPP IPP IP յ թ (եռ յթ յ թ յյթ յյ թ ka 


.000 
1462 
.000 
. 136 
“000 
134 
.000 
-134 
. 000 
126 
.000 
121 


000 


. 126 
2900 
. 131 
. 000 
130 
.000 
1-8 
.000 
125 
.000 
137 
2900 
1595 
551310, 
213 
. 509 
2133 
-500 
2132 
«900 
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299 
1 .000 .000 102.457 -14.000 15.145 14.699 -.016 
5 .018 
2 .000 .250 105.302 -15.000 14.925 14.699 -.013 
.609 -.034 
3 .000 . 500 108.193 -16.000 14.512 14.699 -.015 
. 604 -.092 
4 .000 .750 111.256 -16.000 14.187 14.699 =. 012 
.598 -.152 
5 .000 1.000 112.435 -15.000 13.569 14.698 -.016 
596 -.180 
6 .000 1.250 116.266 -12.000 12.188 14.696 -.044 
. 586 -.288 
7 . 000 1.500 116.796 -5.000 6.156 14.676 -.266 
566 -.522 . 
8 .000 1.750 113.592 12.000 -2.410 1/1: 5341 -.783 
.526 -.970 
9 .000 2.000 105.727 15.000 -18.Կ50 14.573 -1.436 
.483 -1.465 
10 .000 2.250 112.650 72.000 -10.308 14.604 -1.090 
501 -1.258 
Ji .000 2.500 121.443 80.000 . 145 14.651 -.555 
531 -.913 _ 
15 .000 2.750 122.507 76.000 2.894 14.641 -.674 
519 -1.055 
13 .000 3.000 114.758 62.000 -.907 14.573 -1.436 
467 -1.649 ` 
14 .000 3.250 115.506 29.000 -3.502 14.564 -1.543 
. 56 -1.771 
15 .000 3.500 127.261 7.000 2.973 14.639 -. 689 
508 -1.179 ' 
16 .000 3.750 124.883 -6.000 4.305 14.680 24227 
553 -.663 i 
17 .000 4.000 119.450 -11.000 5.371 14.692 -.095 
576 -.409 
18 .000 4.250 114.584 -12.000 5.783 14.694 -.063 
588 -.272 
19 .000 4.500 109.382 -13.000 7.370 14.695 -.059 
598 -.161 
20 .000 &.750 104.913 -13.000 8.771 14.695 -.061 
.605 -.074 
21 .000 5.000 102.151 -13.000 10.026 14.695 -.060 
. 610 -.021 
22 . 000 5.250 99.973 -12.000 10.983 14.694 -.067 
.613 2013 
23 . 000 5.500 99.244 -12.000 11.734 14.694 -.070 
. 6145 .024 
jl .250 .000 99.007 -13.000 16.599 14.689 -.128 
609 -.030 
2 .250 .250 107255 -14.000 15.663 14.694 -.064 
610 -.026 
3 .250 . 500 105.584 -13.000 15.208 14.696 -.049 
605 -.075 
4 .250 .750 109.947 -15.000 14.537 14.695 -.051 
597 -.164 
5 . 250 1.000 111-917 -15.000 13.818 14.696 -.046 
594 -.199 
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